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1.  INTRODUCTION  AND  SUMMARY 

1.1  Introduction 

The  SDI  requirements  for  near  real-time  identification,  discrimination  and  precision 
tracking  of  boosters,  post-boost  vehicles,  and  midcourse  targets  present  a  formidable 
technological  challenge.  As  discussed  below  there  are  many  situations  where  the  performance 
requirements  cannot  be  met  with  conventional  imaging  and  tracking  systems  and  innovative 
techniques  are  needed. 

For  example,  recent  data  on  booster  plume  and  hardbody  signatures  suggest  that  many 
of  the  conventional  passive  and  active  fine  tracking  schemes  are  inadequate  RRI’s  study 
for  STARLAB  indicates  that  the  technical  value  of  near-term  space-based  demonstrations 
must  be  enhanced  through  a  broader-scope  experimental  program  directed  toward  assessing 
the  feasibility  of  precision  tracking  and  pointing  concepts.  Due  to  an  insufficient  data  base  on 
booster  signatures,  the  feasibility  of  crucial  passive-to-active  handover  functions  with 
sufficient  accuracy  is  yet  to  be  demonstrated;  the  hardware  requirements  for  passive  and 
active  trackers  are  yet  to  be  credibly  defined;  and  the  tracking  and  fire  control  algorithms  are 
yet  to  be  developed. 

For  PBV’s  and  midcourse  targets,  there  is  the  well  known  problem  that  conventional 
monostatic  laser  radars  can  be  easily  "blinded”  by  simple  optical  countermeasures  such  as 
comer  cubes.  These  countermeasures  have  forced  the  use  of  bistatic  laser  radar  architectures, 
which  significantly  increase  their  cost  and  complexity.  In  addition,  the  feasibility  of  bistatic 
ladars  for  DEW  applications  has  not  been  evaluated,  particularly  in  terms  of  integration  and 
boresighting  with  the  weapon  platform. 

Typically,  innovative  techniques  are  the  result  of:  1)  Newly  found  properties  of  target 
signatures,  e.g.,  spatio-temporal  fluctuations  of  passive  booster  plume  signatures^’  ,  or 
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radiation  scattering  by  solid-fuel  booster  plumes  in  the  UV,  the  visible,  and  the  near  IR*’4;  2) 
Newly  developed  hardware  components,  e.g.,  coherent  medium-power  lasers  operating  in  the 
UV  or  in  the  visible;  3)  New  algorithms  for  extracting  the  target  information  which  constitute 
the  basis  for  novel  receiver  structures,  e.g.,  laser  correlography,  ,0  imaging  correlography, 
autodyne^  and  image  texture  tracking^,  self-reference  holography  etc. 

RRI  has  made  major  contributions  towards  a  realistic  assessment  of  the  feasibility  of 
active  high-resolution  imaging  for  target  identification  and  precision  tracking  of  interest  to 
SDIO.  Our  research  produced  the  first  quantitative  characterization  of  passive  and  active 
booster  signatures^’4,  which  was  then  used  by  the  SDIO  community  to  predict  the 
performance  of  conventional  fine  tracking  of  boosters.  This  performance  analysis  led  to  the 
conclusion  that  in  numerous  cases  conventional  techniques  are  not  adequate  and  innovative 
approaches  are  necessary. 

Riverside  Research  Institute  proposed  three  novel  techniques  for  fine  tracking  to  meet 
the  unprecedented  performance  requirements:  Autodyne  and  Image  Texture  Tracking  for 
boosters  and  Coherent-Radiation  Imaging  and  Tracking  via  nomination  Coding  (CRITIC)  for 
robust  imaging  and  tracking  of  post-boost  and  midcourse  targets  and  satellites.  These 
techniques  are  described  in  detail  in  Section  3  and  Appendix  l  and  can  be  demonstrated  using 
the  1ST  ground-based  experimental  facilities. 

The  challenge  facing  the  1ST  is  to  develop,  test,  and  select  the  techniques  or  sensor 
systems  for  operational  space  conditions  based  on  experiments  performed  with  the  ground- 
based  sensors.  Thus,  one  will  have  to  develop  the  experimental  field  facilities  of  interest  to 
the  1ST  (the  ISTEF  in  Florida  and  the  Downrange  in  Puerto  Rico)  to  a  state-of-the-art  level 
which  will  allow  them  to  meet  the  following  interrelated  objectives:  1)  to  collect  high- 
resolution  signature  date  on  targets  of  interest  and  2)  to  perform  proof-of-concept  experiments 
with  both  innovative  and  conventional  systems.  To  meet  the  above  objectives  one  has  to 
design  experiments  such  that  their  results  will  be  conclusive  and  scalable  to  operational 
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conditions.  In  Section  2  we  describe  a  novel  ladar  concept,  developed  during  RRI’s 
evaluation  of  1ST  facilities,  that  offers  the  potential  for  real-time  measurements  of 
atmospheric  transmission.  This  technique  will  be  important  for  addressing  scalability  issues. 
Realistic  estimations  of  detection  performance  of  ladars  against  threats  is  dependent  upon 
knowning  atmospheric  transmission  over  the  slant  path  of  the  encounter  geometry. 
Consequently,  atmospheric  extinction  coefficient  is  a  parameter  which  needs  evaluation  to 
support  esdmadons  and  interpretations  of  test  results. 

1.2  Summary 


Collection  of  quantitative  data  on  signatures  of  space  objects  requires  knowledge  of 
atmospheric  transmission  at  the  geographical  site  of  interest.  Many  different  techniques  for 
measuring  atmospheric  transmission  have  been  developed.  In  general,  direct  measurements  of 
transmission  as  a  function  of  altitude  require  an  airborne  transmitter  or  receiver  and  are, 
therefore  difficult  to  cany  out  in  real  time.  Remote  sensing  of  the  atmosphere  by  single-ended 
ladars  avoids  this  difficulty,  but  at  a  price  of  often  unjustified  assumptions  about  extinction 
and  bacKscattering  coefficients.  The  recent  ladar  techniques  which  avoid  such  assumptions 
have  other  types  of  difficulties. 

Review  of  different  techniques  for  transmission  measurements  led  to  the  conclusion 
that  such  measurements  are  either  difficult  to  implement  in  real  time,  or  require  untested 
assumptions  about  the  nature  of  extinction,  scattering,  or  multiple  scattering.  Therefore,  we 
began  the  development  of  a  novel  real-time  multipath  ladar.  The  main  disadvantage  of  this 
multipath  ladar,  common  to  all  remote  sensing  ladars,  is  that  the  signal  may  be  limited  by 
weak  scattering  at  high  altitudes.  On  the  other  hand,  it  has  many  advantages.  It  is  very  simple 
and  except  for  the  steering  mechanism  can  be  built  with  off-the-shelf  hardware.  Measurements 
can  be  made  in  real  time.  Simple  data  processing  is  required  to  obtain  the  values  of 
transmission  coefficient  as  a  function  of  range  without  a  priori  assumptions  about  scattering 
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and  extinction  and  absolute  calibration  is  not  required.  The  multipath  ladar  can  be  tested  by 
varying  the  beam  azimuth  angle  for  a  fixed  zenith  angle.  Radiation  scattered  by  horizontally 
homogeneous  atmosphere  is  independent  of  the  azimuth  angle.  Therefore,  comparisons 
between  the  signals  returned  for  different  azimuth  angles  can  be  used  to  determine  the 
characteristic  length  scales  for  horizontal  inhomogeneities  in  scatterer  properties  at  different 
aldtudes. 

Sizing  of  ground-based  active  systems  for  the  detection  of  space  objects  and  of  the 
proposed  multipath  requires  information  about  the  atmospheric  transmission.  In  the  absence 
of  data  collected  at  a  site  of  interest,  LOWTRAN7  was  used  to  estimate  atmospheric 
transmission  as  a  function  of  range  for  different  wavelengths.  In  addition,  molecular 
contribution  to  the  backscattering  was  calculated  in  the  Rayleigh  limit  and  aerosol 
contribution  using  the  Mie  scattering  for  standard  atmospheric  models.  At  present,  the  validity 
of  these  models  cannot  be  assessed.  A  detailed  discussion  of  this  research  is  provided  in 
Section  2.0. 

In  Section  3  we  evaluate  optical  autodyne  detection  techniques  which  offer  the  potential 
of  improving  the  performance  of  passive  and  conventional  active  sensors  in  detection, 
recognition  and  tracking.  The  autodyne  technique  can  be  applied  against  ground  and 
aerospace  targets  in  a  variety  of  scenarios  as  well  as  for  booster  and  satellite  imaging  and 
tracking.  The  work  presented  here  is  part  of  a  continuing  effort  to  develop  an  analytic 
framework  for  the  detection  and  parameter  estimation  of  autodyne  signals.  The  analysis  is 
applicable  in  a  variety  of  scenarios  involving  rotating  (or  vibrating)  satellites,  or  in  tracking 
the  booster  hardbody-plume  interface.  This  analysis  will  be  useful  in  providing  a  framework 
for  evaluating  other  techniques  based  on:  light  intensity  interferometry  such  as  correlography; 
techniques  applicable  to  diffuse  and  complex  targets  with  strong  specular  components;  and 
will  be  extended  to  the  more  general  case  of  spatio-temporal  light  intensity  interferometry 
where  correlography  and  autodyning  are  limiting  special  cases. 


C49-701 


-4- 


RIVERSIDE  RESEARCH  INSTITUTE 


Optical  autodyne  detection  is  a  direct  detection  method  for  measuring  relative  Doppler 
shifts  using  a  coherent  laser  illuminator.  Various  frequency  components  of  the  received 
optical  signal  are  mixed  together  on  a  square-law  photodetector  to  produce  an  output 
containing  the  spectrum  of  frequencies  related  to  the  relative  velocity  spectrum  of  the  target 
and  background  components.  The  main  advantage  of  autodyne  detection,  as  compared  with 
the  usual  heterodyning  techniques,  is  its  simplicity;  a  local  oscillator  is  not  necessary.  The 
increased  robustness  of  an  autodyne  system  makes  it  attractive  for  many  applications.  Such 
an  active  system  may  consist  of  a  pulsed  laser  transmitter  and  a  photon  bucket  receiver, 
although  an  array  of  autodyne  detectors  in  the  focal  plane  of  the  imaging  sensor  may  also  be 
employed.  One  deficiency  of  autodyne  detection  is  that  it  provides  information  about  the 
autocorrelation  of  functions  proportional  to  the  reflectivity  of  elements  moving  with  a  given 
radial  velocity,  while  the  heterodyne  technique  in  principle  can  provide  direct  information 
about  the  reflectivity  as  a  function  of  radial  velocity.  In  many  scenarios,  however,  a  simple 
autodync  detector  will  provide  a  signature  sufficiently  complementing  the  information 
obtained  by  the  passive  -  active  sensor  suite. 

The  results  of  our  research  are  described  in  detail  in  Section  3  and  summarized  below. 

Fundamental  limits  were  established  for  the  accuracy  of  the  estimation  of  range  and 
velocity  resolution  for  autodyne  detection  of  a  two  point  target  with  a  pulsed  laser  radar 
consisting  of  a  transmitter  collocated  with  a  photon  bucket  receiver.  The  Cramer-Rao  bounds 
(CRB’s)  were  computed  for  both  the  signal  shot  noise  and  background  limited  cases.  From 
these  bounds,  the  number  of  signal  photons  necessary  to  estimate  a  parameter  with  a  desired 
accuracy  may  be  determined. 

The  sensitivity  of  the  bounds  to  different  parameter  values  and  changes  in  the  number 
of  unknown  parameters  was  established.  The  value  of  a  priori  knowledge  in  parameter 
estimation  was  demonstrated  and  related  to  the  number  of  unknown  parameters.  An 
important  issue  of  practical  consequence,  which  remains  to  be  addressed  in  future  work,  is  the 
effect  of  detector  response  time  on  parameter  estimation  of  echo  signals. 
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The  double  point  target  provides  us  with  a  wide  parameter  space.  This  work  would  be 
useful  in  introducing  the  natural  metric  in  parameter  space  for  the  proper  evaluation  of  the 
information  content  in  echo  signals.  The  bound  on  estimator  performance  was  obtained  for  a 
single  laser  pulse.  Evaluating  the  performance  of  a  multipulse  system  is  cenainly  of  great 
interest  and  will  be  addressed  in  future  efforts. 

One  of  the  important  results  obtained  was  that  for  a  large  relative  Doppler  shift  a  few 
tens  of  photons  are  required  to  measure  Doppler  shift  with  moderate  accuracy.  To  put  such 
predictions  on  a  firm  footing,  one  has  to  construct  estimators  with  variances  approaching  the 
CRB.  An  example  applicable  to  problems  in  which  fine  parameter  estimation  is  of  interest, 
e.g.,  alignment  and  calibration  when  the  uncertainty  of  the  parameter  value  is  small,  is 
provided  (e.g,  equation  3.41). 

Another  important  class  of  estimators  useful  for  many  applications  are  the  maximum 
likelihood  estimators  (MLE).  These  estimators  are  applicable  for  large  parameter  variations, 
and  approach  the  CRB  for  large  signal  values.  The  methods  of  implementation  for  the  MLE 
estimators,  together  with  other  practical  considerations,  and  types  of  estimators  are  important 
problems  which  should  be  addressed  in  future  work. 

It  is  expected  that  this  ongoing  effort  will  lead  to  the  development  of  a  new  class  of 
autodyne  sensors  useful  ;n  various  scenarios  of  practical  interest.  It  will  also  help  to  establish 
which  additional  measurements  of  channel,  target,  and  background  signature  are  required,  and 
what  additional  data  related  to  sensor  pattern  selection  are  necessary  to  arrive  at  a  realistic 
assessment  of  the  feasibility  of  advanced  sensors  for  1ST  applications. 
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2.  REAL-TIME  MEASUREMENTS  OF  ATMOSPHERIC  TRANSMISSION 

Knowledge  of  atmospheric  transmission  is  important  in  many  situations.  For  example, 
to  detect  objects  in  the  earth's  atmosphere,  the  wavelength  of  measurement  is  often  chosen  to 
minimize  atmospheric  losses  along  the  line-of-sight  (LOS).  Even  when  detection  is  not  a 
problem  it  may  be  desirable  to  obtain  quantitative  data  on  the  radiation  scattered  or  emitted  by 
the  object  under  study.  In  this  case  radiation  detected  must  be  corrected  for  atmospheric 
losses.  The  values  of  transmission  (or  extinction)  coefficient  as  a  function  of  wavelength 
could  also  be  used  to  determine  concentrations  of  various  molecular  species  in  the 
atmosphere,  and  to  detect  the  presence  of  pollutants  such  as  smoke.  Sizing  of  ground-based 
active  syster.  <  for  the  detection  of  space  objects  also  requires  information  about  transmission 
losses  in  the  atmosphere. 

Not  surprisingly,  many  different  techniques  for  measuring  atmospheric  transmission 
have  been  develop.  These  techniques  are  reviewed  in  Sec.  2.1.  In  general,  direct 
measurements  of  transmission  as  a  function  of  altitude  require  an  airborne  transmitter  or 
receiver  and  are,  therefore,  difficult  to  carry  out  in  real  time.  Remote  sensing  of  the 
atmosphere  by  single-ended  ladars  avoids  this  difficulty,  but  at  a  price  of  often  unjustified 
assumptions  about  backscattering  and  extinction  coefficients.  The  bipath  ladar4’5,  which 
does  not  make  such  assumptions,  requires  an  airborne  ladar  and  again  would  be  difficult  to 
use  in  real  time.  The  multiple-field-of-view  ladar^  requires  its  own  set  of  assumptions 
about  multiple  scattering  in  the  atmosphere,  and  has  not  yet  been  tested. 

Collection  of  quantitative  data  on  signatures  of  space  objects  requires  real-time 
measurements  of  atmospheric  transmission  at  the  geographical  site  of  interest.  In  order  to 
carry  out  such  measurements,  we  began  the  development  of  a  novel  real-time  multipath  ladar 
which  is  described  in  detail  in  Sec.  2.2. 

The  collected  knowledge  of  atmospheric  properties  is  summarized  in  standard 
atmospheric  models®.  Based  on  these  models,  codes  like  LOWTRAN  and  FASCODE  were 
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developed  to  predict  transmission  as  a  function  of  wavelength  and  range  for  a  chosen  LOS; 
some  predictions  of  codes  are  discussed  in  Sec.  2.3.  The  models  attempt  to  describe  "typical" 
atmospheric  conditions,  assumed  to  vary  only  with  altitude,  and  may  be  inappropriate  for  a 
particular  geographic  site.  In  addition,  these  static  models  cannot  account  for  temporal 
variability  of  atmospheric  conditions  which  is  easily  observed  at  low  altitudes.  Finally,  the 
predicted  extinction  coefficient  is  highly  sensitive  to  the  assumed  aerosol  composition^, 
which  is  not  well  known  and  variable  as  well.  Nevertheless,  in  the  absence  of  data  collected  at 
a  geographical  site  of  interest,  the  available  codes  provide  the  means  to  estimate  the  values  of 
the  atmospheric  transmission,  which  are  not  too  far  removed  from  the  realistic  values.  At 
present,  however,  the  validity  of  this  assumption  cannot  be  assessed. 

Sec.  2.4  contains  conclusions  and  recommendations. 

2.1.  Review  of  Techniques  for  Atmospheric  Transmission  Measurements 

Various  techniques  for  measuring  atmospheric  transmission  have  been  reviewed  in 
Refs.  1-3.  They  can  be  divided  into  three  categories:  contrast,  transmission,  and  scattering 
measurements. 

Contrast  measurements  determine  the  apparent  radiance  of  a  known  target  (either 
natural,  such  as  sky,  or  artificial)  and  can,  therefore,  be  used  to  calculate  atmospheric  losses. 
While  very  sensitive  and  portable,  these  "fair-weather"  measurements  can  only  be  made  in 
daylight.  In  addition,  measurements  as  a  function  of  altitude  would  require  airborne  targets. 

Transmission  measurements  can  rely  on  either  natural  or  artificial  sources. 
Measurements  of  starlight  attenuation  determine  only  the  total  optical  depth  of  the  earth 
atmosphere.  Artificial  sources  (lamps  or  lasers)  could  be  used  at  various  distances  from  the 
receivers  to  determine  transmission  as  a  function  of  range.  The  main  difficulty  is  the  proper 
alignment  of  the  transmitter  and  receiver.  However,  these  measurements  are  most  direct  and 
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could  be  used  to  test  other  techniques.  Again,  an  airborne  platform  for  either  the  transmitter 
or  receiver  is  required  for  measurements  at  different  altitudes. 

One  class  of  scattering  measurements,  nephelometry,  relates  the  unknown  extinction 
coefficient  to  the  scattering  coefficients  measured  at  various  scattering  angles.  These 
techniques  work  well  only  if  absorption  is  negligible.  In  addition,  they  require  samples  of  air 
to  perform  measurements,  and  are,  therefore,  difficult  to  use  for  measurements  as  a  function  of 
altitude. 


The  other  class  of  scattering  measurements  uses  single-ended  ladars  in  which  the 
transmitter  and  receiver  are  co- located.  The  power  returned  to  the  receiver  from  radiation 
scattered  at  a  range  R,  from  the  transmitter,  is  given  by  the  standard  ladar  equation, ^ 


P(R)  cc 


"2  Sdx  k  (r) 
0 


(2.1) 


where  P0  is  the  transmitted  power,  0(R)  is  the  volume  backscattering  coefficient,  and  kext(R) 
is  the  extinction  coefficient.  It  is  evident  that  the  measured  signal  P(R)  is  a  function  of  both  0 
and  kcx(;  in  order  to  determine  both  of  these  quantities  from  the  ladar  returns,  it  is  necessary 
to  assume  a  special  relationship  between  the  extinction  kcxf(R)  and  backscattering  0(R). 


It  is  desirable  to  extract  both  0(R)  and  kcxt(R)  from  the  ladar  data  without  assuming  any 
ad  hoc  relationship  between  the  two.  Recently  two  methods  have  been  proposed  for 
achieving  this.  One  is  the  bipath  ladar4’^  and  the  other  is  the  multiple-field-of-view  ladar.^ 


The  bipath  method  employs  two  separate  single-ended  ladar  systems  to  measure  both  the 
backscattering  and  extinction  coefficients  unambiguously.  The  second  ladar  provides  a 
second  independent  equation  in  the  two  variables  of  interest;  this  set  of  two  equations  for  two 
unknowns  can  be  easily  solved.  Hughes  and  Paulson^  have  shown  that  the  backscattering  and 
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extinction  coefficients  obtained  from  a  single  ladar  using  Klett’s  reconstruction  algorithm*® 
were  significantly  different  from  those  obtained  with  the  bipath  method. 

Since  it  is  often  necessary  to  make  measurements  as  a  function  of  altitude,  one  of  the 
ladars  must  be  flown  on  a  platform  (either  a  plane  or  a  balloon).  In  addition,  both  ladars  must 
be  calibrated.  Thus,  the  bipath  method  is  difficult  to  use  for  real-time  measurements  and  a 
single-ended  technique  would  be  preferable. 

One  such  single-ended  method,  the  multiple-field-of-view  ladar,  has  been  suggested  by 
Bissonnette  et  al.^  at  DREV  in  Canada.  This  method  makes  use  of  information  about 
backscattering  and  extinction  which  is  contained  in  mutiply  scattered  radiation  returned  to  the 
receiver.  Although  no  assumption  is  made  about  a  relationship  between  the  backscattering 
and  extinction  coefficients,  Bissonnette  et  al.  do  make  some  assumptions  about  the  multiple 
scattering  contributions  to  the  ladar  signal  which  will  have  to  be  verified.  At  present  this 
system,  built  by  Optech  in  Toronto,  is  being  tested  at  DREV  in  Canada. 

2.2  Novel  Single-Ended  Multipath  Ladar 

As  discussed  in  Sec.  2.1,  ladars  can  sense  atmospheric  properties  in  real  time.  This  great 
advantage  comes  at  a  price  of  often  unjustified  assumptions  about  backscattering  and 
extinction  coefficients.  The  bipath  ladar  which  does  not  make  such  assumptions  requires  an 
airborne  ladar  and  is,  therefore,  difficult  to  use  in  real  time.  The  multiple-ficld-of-vicw  ladar 
has  not  yet  been  tested  and  requires  its  own  set  of  assumptions  about  multiple  scattering  in  the 
atmosphere.  Here  we  describe  a  novel  multipath  ladar  which  can  be  used  in  real  time  and 
which  requires  only  some  weak  assumptions  about  horizontal  inhomogeneities  in  the  earth 
atmosphere. 
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Consider  a  standard  ladar  pointed  along  a  LOS  with  zenith  angle  8  and  azimuth 
angle  8-  This  is  illustrated  schematically  in  Fig.  2.1.  Energy  E(R)  of  backscattered  radiation 
from  the  range  R  detected  by  the  receiver  during  integration  time  is 


T  A  0 

E(X,R)  -  E  — 2_£  $  (X)  {  (R)  0  (X, R) «”2  S*V  1C«xt(X,r) 
O  __Z  0 

2R 


(2.2) 


where  E0  is  the  transmitted  pulse  energy,  A0  is  the  receiver  aperture,  $(X)  is  the  receiver 
transmission  factor,  {(R)  is  the  geometrical  form  factor  which  depends  on  the  transmitter  and 
receiver  geometry,  and  c  is  the  speed  of  light.  Let  Nj  denote  the  number  density  of  scatterers 
(molecules  or  aerosols).  Then  the  volume  scattering  coefficient  is 


0 (X, R)  -  J  H  (R)a  (X^8)  , 
i 


(2.3) 


and  extinction  coefficient  is 


k^lX'Rl  -  2  Ni(R)o*Xt(X) 


(2.4) 


where  ojfX.d5)  is  the  differential  scattering  cross-section  at  an  angle  8s  (usually  for  ladars 
8s  =>  *)  and  OjCxt(X)  is  the  total  extinction  cross-section  which  contains  contributions  from 
both  scattering  and  absorption. 


In  the  coordinate  system  shown  in  Fig.  2. 1 

R  ■  z/co w8 
and 

x  »  ztan 8  . 

Hence,  the  scattering  coefficient  may  be  written  as 

0 (R)  ■  0(*,z)  m  0(ztanfi, z) 
and,  similarly,  extinction  coefficient  becomes 


(2.5) 

(2.6) 

(2.7) 


k  t(R) 
•xt 


k  .  (x,z) 

•xt 


k#xt(zt*n«,  z)  , 


(2.8) 
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Fig.  2.1  Schematic  illustration  of  the  multipath  lidar  for  reai*time  measurements  of 
atmospheric  transmission  and  scattering. 
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where  the  wavelength  dependence  is  implicit.  In  this  coordinate  system  the  ladar  equation 
(2.2)  can  be  written  as 


cr. 


cosd  o  2 


~  COS2 6  5  (X) ( 


cos  0 


0  (stand,  z) 


axp 


{- 


cos  0 


jdz  tand,  z  )J  . 


(2.9) 


Consider  twc  lines-of-sight  with  zenith  angles  d  and  d  +  Ad,  respectively.  Let  Ad  be 
sufficiently  small  so  that 


31nff (x, z) 
3x 


Ax(z) 


«  1 


and 


(2.10) 


31n  ^(x,*) 
3x 


Ax(z) 


«  1  , 


where 


(2.11) 


Ax (z)  2  z (tan(d+Ad)  -  tand]  ,  (2.12) 

Conditions  (2.10)  and  (2.1 1)  indicate  that  the  angular  separation  between  two  beams  must  be 
less  than  the  typical  angular  size  of  horizontal  inhomogeneities  in  the  scatterer  distribution 
and  composition.  These  conditions  are  essential  for  the  proposed  multipath  ladar.  They  are  not 
very  restrictive  since  horizontal  inhomogeneities  are  expected  to  be  most  pronounced  at  lower 
altitudes  where  the  horizontal  separation  between  the  two  beams  is  small.  At  high  altitudes 
the  beam  separation  becomes  large  but  the  atmosphere  is  probably  quite  homogeneous  at 
these  altitudes.  A  technique  for  testing  these  assumptions  is  described  later  in  this  section. 


On  the  other  hand,  let  Ad  be  large  enough  so  that 


_1 _ 1 _ 

cozd  coz(d+Ad) 


>  e. 


(2.13) 


C49-701 


-  14- 


RIVERSIDE  RESEARCH  INSTITUTE 


where  «  is  a  constant  chosen  to  maximize  the  path  length  differences  for  the  two  beams 
subject  to  constraints  (2.10)  and  (2. 11). 

The  values  of  A0  which  satisfy  condition  (2.13)  for  e  =  0.1  are  shown  in  Fig.  2.2.  It  is 
seen  that  the  required  angular  separation  between  the  two  beams  decreases  rapidly  with 
increasing  zenith  angle.  Since,  with  the  exception  of  clouds,  the  characteristic  length  scales 
'  ~  horizontal  inhomogeneities  in  the  atmosphere  are  not  known,  the  multipath  ladar  should 

ate  at  a  zenith  angle  of  about  45°.  Horizontal  separation  between  the  two  beams  is  shown 
in  Fig.  2.3  in  units  of  altitude  as  a  function  of  zenith  angle  9.  It  is  seen  that  for  9  >  45°,  the 
horizontal  separation  becomes  insensitive  to  the  value  of  0;  the  minimum  separation  between 
the  beams  is  ez. 


Assuming  that  conditions  (2.10)  -(2.13)  are  satisfied  we  can  determine  the  extinction 


coefficient  as  a  function  of  altitude  from 


k  (ztan0,z) 
•xt 


i _ i _ 

co »(0+A0)  coaflj  3z  n  *  ' 


(2.14) 


where 


Q<*)  s 


E(z/co b9)  cob  (0+A0)  t (z/coa (0+A0) ) 


E (z/cos ( 0+A0) )  co92(9) 


{ (z/co »9) 


(2.15) 


Since  Q(z)  depends  only  on  the  ratio  of  the  received  signals,  absolute  calibration  is  not 
required  to  obtain  the  extinction  coefficient.  This  is  also  true  for  the  one-way  transmission 
coefficient  which  is  given  by 


T (R*z/co»0)  -  IQ(z)J 


cob (0+A0) _ 

2(co«(0+A0)  -  coa0] 


(2.16) 


Given  the  value  of  the  transmission  coefficient,  the  volume  scattering  coefficient  0  can  be 
obtained  using  the  ladar  equation  (2.2).  Note  that  this  requires  absolute  calibration  of  the 


ladar. 
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The  formulae  (2.14)  and  (2.16)  for  the  extinction  and  transmission  coefficients  were 
derived  under  an  additional  assumption  that  the  scattering  cross-section  depends  on  scattering 
angle  but  not  on  the  incidence  angle.  This  is  true  for  spherical  or  randomly  oriented 
nonspherical  scatterers.  This  assumption  may  break  down  in  the  case  of  nonspherical  particles 
which  may  be  partially  aligned,  for  example  ice  particles  in  cirrus  clouds.  In  such  cases  it  may 
become  necessary  to  carry  out  measurements  at  higher  zenith  angles.  Effects  of  such  possible 
particle  alignment  on  polarization  (or  depolarization)  of  scattered  light  should  be  considered. 

The  multipath  ladar  can  be  tested  by  varying  the  azimuth  angle  4  of  the  beam  for  a  fixed 
value  of  the  zenith  angle  6,  as  shown  in  Fig.  2.1.  Under  the  assumption  of  horizontal 
homogeneity,  the  signal  received  from  a  given  altitude  should  be  independent  of  the  azimuth 
angle  4.  Comparisons  between  received  signals  as  a  function  of  altitude  and  azimuth  angles 
can  be  used  to  determine  the  characteristic  length  scales  of  inhomogeneities  at  different 
altitudes.  The  beam  pointing  parameters  could  then  be  adjusted  in  real  time  so  that  conditions 
(2.10)  -  (2.13)  are  satisfied  during  measurements.  The  multipath  ladar  could  also  measure 
scattering  and  extinction  by  inhomogeneities  such  as  clouds  provided  their  angular  size 
exceeds  the  required  angular  separation  between  the  beams.  In  addition,  by  continuing 
measurements  as  a  function  of  the  zenith  angle,  it  may  be  possible  to  determine  the  two- 
dimensional  distribution  of  scattering  and  extinction.  Temporal  variability  could  be 
determined  on  time  scales  which  will  depend  on  the  laser  pulse  rate  and  the  beam  steering 
time. 

2,3  Calculations  of  Atmospheric  Transmission  and  Backscatterine 

In  order  to  design  the  proposed  multipath  ladar,  it  is  necessary  to  estimate  that  fraction  of 
transmitted  energy  which  will  be  collected  by  the  receiver  as  a  function  of  range.  Clearly,  this 
fraction  depends  strongly  on  the  values  of  the  transmission  and  backscattering  coefficients. 
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All  the  transmission  calculations  discussed  in  this  report  were  obtained  using  a 
modification  of  LOWTRAN7  to  determine  transmission  coefficient  at  a  particular  wavelength 
as  a  function  of  range.  This  is  not  the  intended  use  of  LOWTRAN7,  which  was  designed  to 
calculate  transmission  coefficient  integrated  over  some  spectral  band.  For  losses  in 
propagation  of  laser  beams  it  is  more  appropriate  to  use  FASCODE;  such  calculations  will  be 
carried  out  in  the  future. 

In  the  calculations,  percentage  of  radiation  transmitted  in  a  round  trip  from  the  ground* 
based  transmitter  to  the  monostatic  receiver  has  been  determined  as  a  function  of  target  range 
along  a  line  of  sight  with  the  elevation  angle  of  30°.  Calculations  were  carried  out  at  0.53, 
0.7,  1.06,  3.0,  5.0,  and  10.6  ym.  Mid-Latitude  Summer  Atmospheric  M'  Jel  was  chosen,  as 
most  appropriate  for  the  ISTEF  site. 

Specification  of  the  season  determines  the  tropospheric  and  background  stratospheric 
aerosol  profiles  used  in  calculations.  Volcanic  aerosols  in  the  stratosphere  may  also  contribute 
to  the  transmission  losses  but  have  not  been  included  in  the  calculations  since  it  is  not  clear 
which  model  is  appropriate.  It  should  be  noted,  that  their  effect  would  be  to  reduce  the 
calculated  transmission. 

In  addition  it  was  necessary  to  specify  the  nature  of  the  boundary  layer,  i.e.,  0  to  2  km 
altitude.  For  the  location  of  the  ISTEF  site,  the  NAV  Y  MARITIME  model  was  chosen  as  the 
most  appropriate  one.  In  this  model,  the  air  mass  character  in  the  boundary  layer  can  be  varied 
from  open  ocean  to  strong  continental  influence.  In  order  to  obtain  bounds  on  the  expected 
v  :e  of  transmission,  calculations  were  done  for  both  of  these  extreme  cases. 

The  NAVY  MARITIME  model  also  requires  specification  of  the  meteorological  range 
V  which  is  defined  as 

V  -  3.912 /\xt  f  (2 . 17) 
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where  kext  is  the  total  extinction  coefficient  at  0.55  miti.  To  test  the  sensitivity  of  the 
transmission  coefficient,  calculations  were  done  for  V  =  5  km  (low  visibility)  and  V  =  50  km 
(very  high  visibility). 

The  results  are  shown  in  Figs.  2.4  -  2.9,  in  which  transmission  for  the  continental 
boundary  layer  is  shown  by  solid  curves  while  for  the  open  ocean  boundary  layer  by  dashed 
curves.  It  is  seen  that,  according  to  LOWTRAN7,  essentially  all  of  the  losses  occur  in  the 
first  20  -  30  km  from  the  ground  and  losses  at  3.0  and  5.0  nm  are  significantly  greater  than  at 
other  wavelengths,  due  to  strong  absorption  by  water  and  carbon  dioxide.  For  comparison, 
1976  U.S.  Standard  Atmospheric  Model  has  been  used  to  calculate  transmission  at  10.6  um: 
the  result  is  shown  in  Fig.  2.10.  It  should  be  noted  that  the  calculated  transmission  is  very 
sensitive  to  the  choice  of  atmospheric  model. 

Backscattering  in  the  earth's  atmosphere  is  due  to  molecules  and  aerosols.  Molecular 
contribution  to  the  backscattering  coefficient  was  calculated  using  the  Rayleigh 
approximation  for  the  1976  U.S.  Standard  Atmospheric  Model.  The  aerosol  constituents  of 
the  atmosphere  are  not  known  with  high  accuracy  and,  in  addition,  may  be  highly  variable, 
especially  at  low  altitudes.  In  this  report  we  used  the  aerosoi  models  compiled  in  the 
"Handbook  of  Geophysics  and  the  Space  Environment"  published  by  AFGL.®  The  aerosol 
backscattering  was  calculated  using  the  Mie  theory.  *  * 

The  fraction  of  energy  received  by  the  ladar  was  estimated  by  assuming  the  following 
system  parameters: 


rD  =  1  ms, 

AQ  =  0.1  m^, 

$(X)|(R)  *  0.1. 

Note  that  with  this  integration  time,  the  ladar  would  be  capable  of  range  resolution  of  0.15 
km. 
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Fig.  2.4  Two-way  atmospheric  transmission  at  0.53  um  calculated  with  L0WTRAN7  for 
Mid-Latitude  Summer  Atmosphere  (30°  elevation  angle). 
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Fig.  2.5  Two-way  atmospheric  transmission  at  0.7  ym  calculated  with  LOWTRAN7  for 
Mid-Latitude  Summer  Atmosphere  (30°  elevation  angle). 
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Fig.  2.6  Two-way  atmospheric  transmission  at  1.06  urn  calculated  with  L0WTRAN7  for 
Mid-Latitude  Summer  Atmosphere  (30°  elevation  angle). 
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Fig.  2.7  Two-way  atmoipheric  transmission  at  3.0  urn  calculated  with  LOWTRAN7  for 
Mid-Latitude  Summer  Atmosphere  (30°  elevation  angle).  Note  that  the  vertical 


scale  is  from  0  to  1%. 
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Fig.  2.8  Two-way  atmospheric  transmission  at  5.0  ym  calculated  with  LOWTRAN7  for 
Mid-Latitude  Summer  Atmosphere  (30°  elevation  angle).  Note  that  the  vertical 
scale  is  from  0  to  0.2%. 
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Fig.  2.9  Two-way  atmospheric  transmission  at  10.6  <im  calculated  with  LOWTRAN7  for 
Mid-Latitude  Summer  Atmosphere  (30°  elevation  angle). 


C49-701 


RANGE  (k«) 


RIVERSIDE  RESEARCH  INSTITUTE 


Fig.  2.10  Two-way  atmospheric  transmission  at  10.6  um  calculated  with  L0WTRAN7  for 
1976  U.S.  Standard  Atmosphere  (30°  elevation  angle). 
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The  results  are  shown  in  Figs.  2.1 1  and  2.12  for  wavelengths  of  0.53  and  1.06  *tm, 
respectively.  In  these  figures,  the  energy  fraction  is  given  as  a  function  of  altitude  rather  than 
range,  which  depends  on  the  zenith  angle  8.  The  weak  signal  from  higher  altitudes  may 
require  longer  integration  times,  larger  receiver  apertures,  and/or  more  powerful  laser 
transmitter. 

2.4  Conclusions  and  Recommendations 

Collection  of  quantitative  data  on  signatures  of  space  objects  requires  real-time 
measurements  of  atmospheric  transmission  at  the  geographical  site  of  interest.  Review  of  the 
different  techniques  for  transmission  measurements  led  to  the  conclusion  that  either  such 
measurements  are  difficult  to  implement  in  real  time,  or  that  they  require  untested 
assumptions  about  the  nature  of  extinction,  scattering,  or  multiple  scattering.  Therefore,  we 
began  the  development  of  a  novel  real-time  multipath  ladar. 

The  main  disadvantage  of  the  multipath  ladar  is  common  to  all  remote  sensing 
ladars;  i.  e.,  it  relies  on  backscattering  to  provide  the  signal  and  may  thus  be  limited  by  weak 
scattering  at  high  altitudes.  Therefore,  averaging  over  several  pulses  may  be  required  to 
improve  the  SNR. 

One  the  other  hand,  the  multipath  ladar  has  many  advantages.  First  of  all  it  is  very 
simple:  except  for  the  steering  mechanism,  it  can  be  built  with  off-the-shelf  hardware. 
Measurements  can  be  made  in  real  time.  Simple  data  processing  is  needed  to  obtain  the 
values  of  transmission  coefficient  as  a  function  of  range  or  extinction  coefficient  as  a  function 
of  altitude  and  zenith  angle  without  prior  assumptions  about  scattering  and  extinction. 
Determination  of  the  backscattering  coefficient,  however,  requires  absolute  calibration  of  the 
system. 
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Fraction  of  transmitted  energy  collected  by  the  lidar  receiver  as  a  function  of 
altitude  for  different  zenith  angles  and  X  *  0.53  am. 
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Fig.  2.12  Fraction  of  transmitted  energy  collected  by  the  lidar  receiver  as  a  function  of 
altitude  for  different  zenith  angles  and  X  *  1.06  «m. 
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The  multipath  ladar  can  be  used  to  determine  the  characteristic  length  scales  for 
horizontal  inhomogeneities  in  scatterer  properties  as  a  function  of  altitude  by  varying  the 
beam  azimuth  angle  for  a  fixed  zenith  angle.  Extinction  and  scattering  by  inhomogeneities  of 
angular  size  greater  than  the  angular  separation  between  two  beams  could  also  be  measured. 
Such  inhomogeneities  may  include  clouds  and  aircraft  contrails. 

If  a  dye  laser  is  used  in  this  ladar,  the  wavelength  dependence  of  scattering  and 
extinction  coefficients  could  be  determined.  This  is  very  important  for  measuring  molecular 
absorption  and,  therefore,  concentration  of  absorbing  species,  and  also  for  determining  the 
aerosol  size  distribution. 

Several  issues  pertaining  to  the  multipath  ladar  still  have  to  be  addressed.  One  of 
them,  common  to  all  ladars,  is  the  effect  of  clear-air  turbulence  which  may  induce  beam 
wander  or  spreading.  12-14  Narrow  (compared  to  the  turbulence  length  scale)  beams  may 
wander,  leading  to  uncertainties  in  the  zenith  angle,  while  wide  beams  may  spread,  leading  to 
uncertainty  in  the  geometric  form  factor  for  the  ladar.  The  choice  of  the  optimal  beam  width 
must  be  investigated.  Averaging  over  pulses  may  be  useful  in  separating  the  effect  of  clear-air 
turbulence  and  should  also  be  analyzed. 

The  other  issue  is  the  possible  alignment  of  nonspherical  scatterers,  e.  g.,  ice 
particles  in  cirrus  clouds.  Possible  advantages  of  polarization  sensitive  measurements  should 
be  investigated. 

The  proposed  single-ended  multipath  ladar  can  be  tested  by  using  a  transmissometer 
with  an  airborne  receiver  (or  transmitter).  Transmissometer  measurements  are  very  important 
not  only  for  determining  the  overall  accuracy  of  the  proposed  system  but  also  for  assessing  its 
performance  under  different  meteorological  conditions. 
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3.0  AUTODYNE  DETECTION 


1  ^ 

RRI  proposed  to  use  the  optical  autodyne  detection  *•"  system  to  improve  the 
performance  of  passive  and  conventional  active  sensors  in  detection,  recognition  and  tracking. 
The  proposed  technique  can  be  applied  against  ground  and  aerospace  targets  in  a  variety  of 
scenarios  as  well  as  for  booster  and  satellite  imaging  and  tracking.  The  work  presented  here  is 
a  pan  of  the  continuing  effon  to  develop  an  analytic  framework  for  the  detection  and 
parameter  estimation  of  autodyne  signals,  an  analysis  applicable  in  a  variety  of  scenarios  with 
targets  such  as  rotating  (or  vibrating)  satellites,  or  tracking  hardbody-plume  booster  interface. 
This  analysis  will  be  useful  in  evaluating  other  techniques  based  on  light  intensity 
interferometry  such  as  correlography,  for  diffuse  and  complex  targets  with  strong  specular 
components  and  will  be  extended  to  spatio-temporal  light  intensity  interferometry  where 
correlography  and  autodyning  are  limiting  special  cases. 

When  an  incoherent  laser  illuminator  is  replaced  by  a  coherent  one,  active  images  are 
degraded  by  laser  speckle.  It  is  possible,  however,  to  obtain  Doppler  measurements  providing 
information  about  the  target  radial  velocity  and  possibly  its  vibration.  This  signature  permits 
the  detection  of  moving  or  vibrating  targets  in  the  presence  of  fixed  or  static  targets,  even 
when  the  signature  from  the  fixed  target  is  ord:  of  magnitude  greater. 

One  can  employ  heterodyne  reception  in  the  image  plane  when  the  detector  array  is 
illuminated  by  a  local  oscillator  laser  and  outputs  are  properly  processed.  This  provides 
velocity  information  in  each  resolution  cell  of  the  active  imaging  sensor  and  reduces  the 
contribution  from  background  (clutter)  and  detector  noise.  Measurement  of  the  velocity 
distribution  within  the  illuminating  beam  can  also  be  obtained  without  resorting  to  the 
imaging  system,  by  employing  a  single  detector  properly  matched  to  the  size  of  speckle 
produced  by  reflections  of  the  illuminating  beam  and  using  heterodyne  detection.  The 
heterodyne  deter^on  technique  requires  very  precise  frequency  matching  of  illuminator  and 
local  oscillator  lasers  and  precise  alignment  of  the  signal  and  local  oscillator  beams.  This  may 
be  quite  difficult  to  achieve  and  maintain,  especially  for  shorter  wavelengths. 
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RRI  proposed  to  replace  the  sophisticated  and  sensitive  heterodyne  method  by  a 
more  robust  autodyne  technique.  Optical  autodyne  detection  is  a  direct  detection  method  for 
measuring  the  relative  Doppler  frequency  shifts  using  a  coherent  laser  illuminator.  Various 
frequency  components  of  the  received  optical  signal  are  mixed  together  on  a  square-law 
photodetector  to  produce  an  output  containing  the  spectrum  of  frequencies  related  to  the 
relative  velocity  spectrum  of  the  target  and  background  components.  The  main  advantage  of 
autodyne  detection,  as  compared  with  the  usual  heterodyning  techniques,  is  its  simplicity;  a 
local  oscillator  is  not  necessary.  The  increased  robustness  of  an  autodyne  system  makes  it 
attractive  for  many  applications.  Such  an  active  system  may  consist  of  a  pulsed  laser 
transmitter  and  a  photon  bucket  receiver,  although  an  array  of  autodyne  detectors  in  the  focal 
plane  of  the  imaging  sensor  may  also  be  employed.  One  deficiency  of  autodyne  detection  is 
that  it  provides  information  about  autocorrelation  of  functions  proportional  to  the  reflectivity 
of  elements  moving  with  a  given  radial  velocity,  while  heterodyne  technique  in  principle  can 
provide  direct  information  about  the  reflectivity  as  a  function  of  radial  velocity.  In  many 
scenarios,  however,  a  simple  autodyne  detector  will  provide  a  signature  sufficiently 
complementing  the  information  obtained  by  the  passive  -  active  sensor  suite. 

A  limited  experimental  proof-of-principle  laboratory  demonstration  has  been  already 
performed^. 

In  the  next  funded  program  phase,  it  is  recommended  that  RRI  design  more  complete 
and  realistic  experiments  illustrating  the  capabilities  of  the  autodyne  technique  for  pulsed 
lasers,  and  define  hardware  requirements  for  the  autodyne  system  in  various  scenarios  of 
interest.  Different  algorithms  capable  of  extracting  relevant  information  for  complex  targets, 
containing  a  mixture  of  diffuse  and  specular  reflectivity  components,  could  be  proposed  and 
sensitivity  studies  performed  to  evaluate  the  applicability  of  such  systems  to  the  ONR 
mission. 


The  program  presented  by  RRI  will  lead  to  the  development  of  a  new  class  of 
autodyne  sensors  useful  in  various  scenarios  of  practical  interest.  It  will  also  help  to  establish 
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which  additional  measurements  of  target,  and  background  signature,  and  what  additional  data 
related  to  sensor  pattern  selection  are  necessary  to  arrive  at  a  realistic  assessment  of  the 
feasibility  of  the  proposed  technique. 

3.1  Backeround,  Technical  Approach,  and  Anticipated  Benefits 

3.1.1  Background 

In  many  applications,  an  improved  performance  of  optical  passive  and  active  sensors 
against  ground,  airspace  and  space  targets  is  required.  This  enhancement  can  be  achieved  by 
the  addition  of  a  coherent  laser  illuminator.  If  an  active  imaging  receiver  is  employed, 
speckled  images  of  the  target  can  be  obtained,  and  with  appropriate  wavefront  modulation  one 
obtains  range  information.  Employment  of  the  autodyne  array  for  an  active  imager  will 
provide  information  about  the  relative  velocity  distribution  in  a  resolution  cell  of  the  active 
imager.  A  photon  bucket  receiver  with  autodyne  detector  will  provide  information  about  the 
relative  velocity  distribution  in  t^illuminated  region.  An  array  of  properly  distributed 
photon  bucket  receivers  will  not  om>  improve  the  appropriately  defined  signal*to-noise  ratio, 
but  may  give  some  information  about  the  spatial  distribution  of  the  target  and  background 
components  and  their  velocities  (imaging  and  Doppler  measurement  with  spatio-temporal 
intensity  interferometry).  While  coherent  (heterodyne)  detection  can  provide  unambiguous 
information  about  absolute  velocity  distribution,  the  autodyne  mode  of  detection  can  offer  a 
practical  scheme  of  sufficiently  robust  and  inexpensive  sensor  for  many  applications. 

A  meaningful  evaluation  of  an  active  system  requires  a  good  understanding  of  the 
role  it  has  to  fulfill  in  various  missions  and  scenarios.  It  could  vary  from  detection  of  a 
moving  target  within  a  set  of  possible  targets  selected  by  passive  sensors  in  a  cluttered 
background  to  the  determination  of  the  type  of  target  from  vibration  spectrum  or  the  detection 
and  tracking  of  a  selected  feature  in  a  complex  target.  The  capability  of  providing  rapid 
illuminator  beam  steering  may  be  required;  weather  constraints  may  be  an  important  factor  to 
be  considered.  Important  elements  defining  optical  system  performance  include  illuminator 
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power,  beam  wavefront  quality,  pulse  repetition  rate,  constraints  and  effects  of  the  selection  of 
illuminator  and  receiver  platform,  the  properties  of  bidirectional  propagation  channel,  target 
and  background  effects,  detector  characteristics,  resolution  and  s  nsitivity  of  spectrum 
analyzers  or  digitizers.  RR1  is  aware  of  these  effects  and  has  the  c  pability  to  bound  and 
quantify  their  contribution  to  the  overall  system  performance.  Understanding  these  issues  and 
their  impact  on  the  performance  of  heterodyne  and  autodyne  systems,  together  with  the 
knowledge  of  the  state  of  the  an  is  essential  to  assess  the  potential  and  the  deficiencies  of 
these  two  techniques. 

3.1.2  Technical  Approach 

Let  us  assume  that  both  receiver  and  transmitter  are  located  on  the  ground.  The 
performance  of  a  highly  coherent  laser  transmitter  illuminating  the  target  with  known 
intensity  distribution  may  be  affected  by  platform  vibration  and  optical  path  phase  and 
amplitude  distortions.  Platform  vibrations  are  important  for  a  sufficiently  narrow  beam,  when 
the  induced  beam  jitter  becomes  comparable  with  the  beam  size.  Part  of  the  optical 
degradation  is  due  to  the  index  of  refraction  fluctuations  caused  by  clear  weather  turbulence^ 
of  the  atmosphere.  The  size  of  these  fluctuations  depends  on  weather  conditions,  transmitter 
and  receiver  location,  transmitter  line  of  sight  (LOS),  and  for  fixed  aperture  size  their 
contribution  increases  as  the  wavelength  becomes  shorter. 

The  presence  of  the  turbulence  effects  may  produce  additional  beam  jitter  or  more 
serious  degradation  of  illuminating  beam  for  larger  apertures  or  ranges. 

It  is  useful  to  divide  the  effects  of  turbulent  atmosphere  into  two  pans  corresponding 
to  segments  of  a  bidirectional  earth  -  space  -  eanh  propagation  channel.  The  highly  coherent 
earth  based  laser  transmitter  illuminates  target  in  space  with  a  beam  of  optically  flat 
distribution  or  other  known  intensity  distribution  (uplink).  Diffuse  and  specular  reflection  of 
the  laser  illuminator  is  detected  by  the  ground  based  receiver  (downlink). 
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A  rudimentary  analysis  reveals  that  turbulence  -  induced  beam  spreading  will 
increase  the  diffraction  limited  far  field  beamwidth  of  the  laser  transmitter,  and  angle-of- 
arrival  fluctuations  may  limit  the  power  collection  capability  of  the  diffraction  limited  receiver 
(if  it  is  sufficiently  large).  While  atmospheric  beam  spreading  may  degrade  system 
performance  operating  in  low  light  conditions  (e.g.,  at  large  ranges),  both  in  correlography  in 
autodyne  detection  (with  nonimaging  receiver),  the  use  of  diffraction  limited  receivers  is  not 
required  and  atmospheric  limitations  on  the  size  of  PSF  are  not  especially  significant  for  a 
receiver  that  is  essentially  an  aperture  integrator  (photon  bucket).  Both  uplink  transmission 
and  downlink  reception  are  subject  to  scintillation  fading,  and  those  effects  may  significantly 
affect  the  performance  of  various  techniques  employing  analysis  of  spatio-temporal  behaviour 
of  speckle  patterns  to  obtain  the  desired  target  characteristic.  In  general,  both  uplink  and 
downlink  effects  will  be  present  simultaneously.  It  is  instructive,  however,  to  consider 
degradations  produced  in  a  simplified  case  when  only  one  of  these  two  effects  is  present.  In 
particular,  compensating  optics  can  be  used  to  combat  uplink  fading,  while  on  the  downlink 
this  can  be  accomplished  by  spatial  averaging  of  atmospheric  speckle  with  sufficiently  large 
receivers  (if  target  generated  speckle  is  sufficiently  large),  or  by  increasing  the  number  of 
receivers. 

If  we  assume  that  the  target  is  far  above  the  atmosphere,  then  all  turbulence  in  the 
propagation  path  is  located  near  the  receiver  and  transmitter  plane.  Because  of  this 
concentration  of  turbulence  near  receiver  and  transmitter,  the  earth-space  channel  behaves  as 
though  the  vacuum  channel  was  dominating  through  mos:  of  the  propagation  path  from 
transmitter  (receiver)  to  target,  with  turbulence  creating  a  random  phase  and  amplitude  screen 
in  the  transmitter  (receiver)  plane.  This  simple  modeling  which  neglects  the  thickness  of 
atmospheric  layer  may  occasionally,  depending  on  applications,  be  inadequate  or  even  lead  to 
contradictions.  Defining  the  conditions  when  simple  modeling  of  the  propagation  channel 
with  one  of  a  few  properly  positioned  random  screens  is  adequate  for  a  given  application 
would  greatly  contribute  to  better  understanding  of  the  requirements  for  the  different 
techniques  of  spatio-temporal  speckle  processing  to  provide  satisfactory  estimates  of  target 
parameters. 
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It  can  be  easily  shown  that  downlink  effects  are  equivalent  to  multiplicative  random 
modulation  intensity  distributions  of  target  generated  speckle  patterns  in  the  receiver  plane 
(provided  that  the  target  is  smaller  than  the  isoplanatic  patch  size).  Similarly,  uplink  effects 
are  equivalent  to  random  modulation  of  illuminator  beam  wavefront  which  multiplies  target 
reflectivity  amplitude.  The  result  of  uplink  and  downlink  turbulence  effects  can  be  presented 
as  generating  different  type  of  bias  of  the  estimates  of  various  target  characteristics. 

It  is  worth  noting  that  the  resolution  of  passive  and  active  sensors  will  be  similarly 
affected  by  turbulence.  If  the  angle-of-arrival  fluctuations  cause  a  degradation  of  active 
images,  this  will  affect  the  autodyne  as  well  as  the  heterodyne  detection  techniques.  If  a  non¬ 
imaging  receiver  that  is  essentially  an  aperture  integrator  (photon  bucket)  is  employed,  this 
atmospheric  degradation  is  not  important  for  the  autodyne  technique,  but  it  may  degrade  the 
efficiency  of  heterodyne  detection  if  the  correladon  scale  of  the  angle-of-arrival  fluctuations  is 
comparable  to  the  detector  size.  If  non-imaging  arrays  of  detectors  are  employed,  the 
presence  of  the  channel  induced  degradation  further  diminish  the  attractiveness  of  heterodyne 
schemes.  Since  both  parts  of  the  propagation  channel  are  subject  to  scintillation  fading 
effects,  they  may  influence  the  performance  of  autodyne  as  well  as  heterodyne  detection 
schemes.  In  some  scenarios,  non-imaging  arrays  of  autodyne  detectors  could  produce  outputs 
which,  when  properly  processed,  would  allow  resolution  exceeding  the  limitations  imposed 
by  the  turbulence  effects  on  imaging  sensors. 

All  this  indicates  that  the  presence  of  atmospheric  turbulence  induced  degradation 
either  worsens  both  systems  in  a  similar  fashion  or  makes  the  autodyne  scheme  more 
attractive.  Other  changes  in  the  signal,  due  to  absorption  by  atmospheric  gases  or  scattering 
by  molecules  or  aerosols  which  could  be  of  great  importance  in  some  scenarios,  will  affect 
both  techniques  in  a  similar  fashion  and  could  be  the  decisive  factors  in  wavelength  selection. 

The  laser  radar  signatures  of  the  different  targets  and  their  proper  characterization  in 
various  operational  scenarios  is  an  important  element  of  performance  evaluation  of  an  active 
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sensor.  Many  targets  of  interest  represent  geometrically  complex  objects.  For  instance,  some 
satellites  should  be  modeled  by  several  flat  and  curved  surfaces  and  edges.  This  may  result  in 
rather  complicated  signatures,  strongly  dependent  on  the  relative  configuration  of  the 
transmitter,  receiver,  and  target.  The  relative  contribution  to  the  observed  signatures  of  the 
specular  component  of  the  target  reflectance  due  to  edges  and  flat  surfaces  normal  to  the 
transceiver  vs  the  diffuse  component  of  the  target  reflectance,  their  stability  and  classification 
potential  will  affect  the  choice  and  mode  of  operation  of  an  active  sensor. 


Let  us  assume  that  the  target  is  illuminated  with  a  beam  which  is  larger  than  the 
target  dimensions,  and  that  longitudinal  coherence  of  the  transmitted  illumination  pulse  is 
determined  by  the  finite  duration  of  the  pulse.  We  may  be  interested  in  target  vibration 
spectrum  where  the  target  may  represent  a  helicopter  in  the  vicinity  of  clouds. 


Let  us  assume  that  Ag  is  an  effective  illuminator  beam  area  in  the  target  plane.  Then 
a  properly  matched  receiver  at  distance  R  should  have  the  dimensions 

Ag  -  (x2/  A^R2.  (3-D 

The  number  of  photons/sec  incident  on  the  receiver  is  given  by 


where  Npj  is  the  number  of  photons/sec  transmitted,  X-  is  laser  wavelength,  and  dO  is  the 
solid  angle  given  by  Ag/R^.  We  also  assumed  that  the  illuminated  spot  can  be  described  by  a 
Lambertian  scatterer  with  effective  reflectivity  <  p.  Therefore, 

"fr-  (hHK'SlK* 


■  |s.l024j [*p/*j (*3/*b|Pt*  (3.3) 

where  Fj  is  the  power  in  watts. 
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The  Doppler  shift  of  backscattered  radiation  is  given  by: 

Av  -  |2v/cJvo»  2v^/X  (3.4) 

where  vQ  is  the  light  frequency,  and  vT  is  the  radial  velocity  of  the  scatterer.  The  smallest 
temporal  scale  is  related  to  the  largest  relative  velocity  Avrmax  present  among  scatterers  (the 
so-called  temporal  size  of  speckle),  and  is  given  by 


T 

•P 


1/Av 


max 


X/2Av 


max. 


(3.5) 


This  important  parameter  defines  the  required  detector  response.  The  number  of  photons 
registered  in  this  characteristic  time  is  given  by 


*  ■  (,/’](x4/ABllVAVrm«*]|5,1°24H1/’'*tt'"*t*t)  (3'6) 

In  order  to  achieve  autodyne  detection  which  is  not  shot  noise  limited  (and  assuming 
that  detectors  are  shot  noise  limited),  a  minimum  N  of  5  ♦  10  photons  are  required  per 
temporal  speckle5;  therefore,  to  operate  in  this  regime,  the  expression  for  the  power 
requirement  is  given  by: 


PT[W] 


2w 

< 

P 


[matar] 


(3.7) 


where  tjt0t  describes  the  total  transmitter  and  receiver  efficiency,  and  includes  possible 
atmospheric  losses. 


The  amount  of  energy  per  pulse  is  determined  by  the  velocity  resolution  requirement 
and  defined  by  pulse  duration.  This  second  characteristic  time  scale  is  given  by 


T 

P 


■  X/Av 


min. 


and  leads  to  the  following  energy  requirement: 


(3.8) 
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*U3  -  ^7 
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max 


X3Av  .  T)  _ 
min  ToT 


f.  ,rt241 
5  •  10 


[mater] 


(3.9) 


This  amount  of  energy  (for  N  of  the  order  of  5  +  10)  would  allow  for  "semiclassical" 
operation.  Note  that  in  addition  to  the  expected  dependence  on  illuminated  area,  receiver  and 
transmitter  efficiency,  it  is  proportional  to  the  inverse  of  the  cube  of  the  wavelength  and  the 
number  of  velocity  resolution  cells  defined  by 


N 


ras 


Av  /Av  , 
max  min 


(3.10) 


Increasing  laser  energy  beyond  that  required  for  semiclassical  operation  will  not  bring  about 
any  significant  improvement  (for  shot  noise  limited  detectors)  and  fixed  pulse  deviation.  In 
semiclassical  operation,  after  proper  signal  processing,  we  can  obtain  autocorrelation  of 
reflected  intensity  disbribution  as  a  function  of  radial  velocity.  Each  "resolution"  cell  in  this 
image  will  be  comipted  by  noise.  Typically,  SNR  defined  as  signal  divided  by  the  square 
root  of  variance  of  the  signal  can  be  of  the  order  of  1/2  (results  depend  on  the  signal 
processing  scheme  and  on  target  signatures). 


Consider  the  following  numerical  example: 


P 

N 


■  0.1, 


«  10, 


"tot  ■  °-05' 

(^max^^min  j 


100  m2, 
100. 


X  -  10 


(3.11) 


This  leads  to  Ej  *  25J.  In  the  visible  range,  for  X  *  .5  •  lO'^m,  Ej  *  200J  and  for  long 
wavelength  X  »  10*  10'^m,  Ej  *  (1/40)J.  This  rather  considerable  energy  requirement  for 
short  wavelength  cannot  in  general  be  easily  decreased  by  any  significant  amount  without 
paying  a  price  somewhere  else.  For  instance,  a  decrease  of  the  energy  by  a  factor  of  two  can 
be  compensated  to  produce  the  same  SNR  by  using  two  autodyne  detectors  or  two  (so  that  the 
number  of  photons  per  temporal  speckle  N  is  reduced  from  10  to  5)  pulses;  further  decrease  of 
the  energy  transmitted  may  require  a  disporportional  increase  of  the  number  of  autodyne 
receivers.  For  lowlight  conditions,  this  number  would  be  proportional  to  the  square  of  the 
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inverse  of  transmitted  energy.  For  longer  wavelength  the  energy  requirements  as  listed  do  not 
represent  a  significant  constraint,  but  the  formulas  are  too  optimistic,  since  infrared  detectors 
are  not  shot  noise  limited.  One  can,  however,  easily  replace  the  formula  presented  above  by 
interpolating  the  expression  which  takes  into  account  the  detector  noise.  Let  us  note  that 
energy  requirements  for  heterodyne  reception  are  similar,  except  that  detector  noise  is  usually 
negligible. 


When  tracking  booster-plume  interface  with  a  narrow  beam  with  effective  area  of 
4m^  we  obtain  much  more  moderate  energy  requirements  in  the  visible  range  of  Ej  =  8J. 


It  is  worthwhile  to  compare  the  energy  requirement  of  non-imaging  autodyne 
detection  and  that  of  an  imaging  sensor  (in  non-autodyne  mode,  to  assure  a  performance 
which  is  not  shot  noise  limited).  In  this  case  the  number  of  photons  incident  on  the  imaging 
receiver  is  given  again  by: 


N. 


[! 


«  /» 
P 


dO 


N. 


[(*/’)  (vR*]h 


(3.12) 


where  Ar  is  the  area  of  imaging  receiver  and  Nj  is  the  number  of  photons  transmitted.  Since 
the  area  of  resolution  cell  in  the  target  plane  is  given  by  (X^/A^)R^,  the  number  of  photons 
from  the  resolution  cell  of  the  illuminated  spot  is  given  by 

n»  ■  [(V'lK'slK 

■  js-io24) [«  /») (x3/AB]ST[l/Joul«.m«t0r]  (3.13) 

/\ 

where  Nj  is  pulse  energ>  in  Joules.  In  order  to  achieve  an  image  quality  which  is  not  shot 
noise  limited  (and  assuming  coherent  illumir  nd  shot  noise  limited  detectors),  we  need 

an  Njr  of  the  order  of  5  per  resolution  cen.  jre,  to  operate  in  this  regime,  we  obtain 

the  following  energy/pulse  requirement: 
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and  t?T0T  again  describes  the  transmitter  and  receiver  efficiency  together  with  atmospheric 
channel  losses. 

This  expression  is  similar  to  the  one  obtained  for  autodyning,  except  that  it  does  not 
involve  a  multiplier  equal  to  the  number  of  velocity  resolution  cells.  Note  that  a  factor  of  two 
is  also  missing  and  our  estimate  of  Njr  necessary  to  attain  the  classical  regime  can  also  be 
lower  by  a  factor  of  two,  one  has  to  be  careful,  however,  with  such  order-of-magnitude 
estimates. 

Thus,  accepting  the  numbers  from  the  previous  numerical  example,  we  have  an 
energy  requirement  which  is  more  than  a  hundred  times  lower  than  in  autodyning.  In 
addition,  for  the  imaging  sensor  we  can  trade  directly  the  laser  energy/pulse  for  the  number  of 
pulses  or  the  number  of  receivers,  which,  as  mentioned  above,  is  possible  only  for  a  limited 
energy  range  in  autodyning.  Another  condition  necessary  for  achieving  a  performance  which 
is  not  shot  noise  limited  is  the  requirement  for  the  minimum  power.  Decreasing  laser  power 
and  increasing  the  pulse  duration  with  fixed  energy  increases  the  resolution  of  the  autodyne 
technique.  For  an  imaging  system  and  low  light  level  conditions,  the  resolution  can  be  traded 
directly  for  SNR.  Similar  direct  tradeoffs  for  autodyning  are  possible  only  for  a  limited  range 
of  power,  without  causing  SNR  degradation.  If,  however,  energy  beyond  that  required  for 
semiclassical  operation  is  available,  increasing  pulse  duration  is  often  beneficial.  The  energy 
requirement  for  an  autodyne  array  of  an  imaging  system  is  similar  to  that  of  a  nonimaging 
system.  It  is  important  to  note  that  employment  of  an  imaging  sensor  with  too  high  resolution 
may  rule  out  the  appearance  of  an  autodyne  signal  if  the  target  signal  components  with 
different  velocities  are  separated  into  different  resolution  cells. 

This  general  discussion  needs  refinement  and  should  be  more  strongly  connected  to 
various  scenarios  and  detected  target  signatures.  For  instance,  when  the  velocity  spectrum 
contains  only  two  points  as  in  the  case  of  aircraft  moving  in  a  clouded  background,  detection 
and  estimation  of  the  signal  may  require  lower  laser  power  than  in  the  case  of  the  general 
velocity  distribution.  This  should  be  included  in  the  discussion  of  the  value  of  a  priori 
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knowledge  in  estimating  the  parameters  of  interest.  A  major  section  of  this  report  is  devoted 
to  this  problem.  Another  important  case  corresponds  to  the  target  with  specular  reflection.  In 
some  cases  glint  reflection  would  serve  as  a  local  oscillator  for  the  reflected  diffuse 
component  and  effectively  convert  the  autodyne  system  into  heterodyne  detection**. 
Improvement  of  SNR  by  the  employment  of  several  autodyne  detectors  was  mentioned  above 
in  the  context  of  shot-noise  reduction.  In  many  cases  using  sparse  arrays  of  autodyne 
detectors  would  lead  to  noise  reduction  due  to  fading.  To  introduce  proper  spatio-temporal 
sampling  redundancy,  helpful  in  improvement  of  SNR  (if  possible),  it  is  necessary  to  have  a 
good  understanding  of  target  signatures,  and  in  particular  the  target  coherence  properties. 

Another  problem  not  discussed  here  is  the  use  of  proper  wavefront  modulation  to 
enhance  the  range  estimation,  and  also  to  obtain  information  about  the  target  depth.  While 
range  information  can  be  easily  extracted  for  relatively  short  pulses,  with  simple  amplitude 
modulation  more  information  about  the  target  in  the  line  of  sight  direction  can  be  obtained 
using  autodyne  detection  by  appropriate  phase  modulation  similar  to  the  one  employed  with 
conventional  radar.  We  will  briefly  touch  upon  this  problem  in  the  next  section. 

3.1.3  Anticipated  Benefits 

The  discussion  of  issues  outlined  above  could  be  carried  out  in  more  detail  for 
various  scenarios  of  interest  in  the  future  to  establish  when  the  autodyne  system  can  enhance 
the  performance  of  a  passive  and  conventional  active  sensor  suite.  The  development  of  an 
analytical  framework  would  help  to  understand  the  potential  of  the  autodyne  system.  This 
would  also  help  to  assess  the  usefulness  of  the  autodyne  photon  bucket  detector  array  in 
obtaining  both  target  images4'^  and  velocity  information  for  various  scenarios.  Such  an 
analysis  will  help  to  select  proper  sensor  architecture  (imaging  or  nonimaging)  for  limited 
demonstration  and  short-  and  long-term  applications.  The  results  obtained  would  be  useful  to 
guide  technology  development  in  the  area  of  lasers,  spectrum  analyzers  and  digitizers. 
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3.2  Fundamental  Limits  on  Ranee  and  Velocity  Resolution  for  Autodvne  Detection 

of  Two  Point  Target 

In  this  section  we  will  consider  limits  on  range  and  velocity  resolution  for  autodyne 
detection  of  a  two  point  target  with  a  pulsed  laser  radar  consisting  of  a  transmitter  collocated 
with  a  photon  bucket  receiver.  In  this  type  of  system  the  target  is  illuminated  with  a  beam 
which  is  typically  larger  than  the  target  dimension  and  the  received  signal  is  registered  by  a 
square-law  photodetector.  The  longitudinal  coherence  of  the  transmitted  illumination  pulse  is 
determined  by  the  finite  duration  of  the  pulse.  In  our  analysis  we  assumed  that  the  amplitude 
of  the  transmitted  pulse  is  Gaussian.  The  considered  model  includes  linear  modulation  of  the 
phase  of  the  transmitted  pulse. 

We  will  discuss  the  estimation  of  ranges,  relative  velocity,  and  amplitudes  of  a  two 
point  target  when  laser  backscattered  radiation  is  collected  with  a  photon  bucket  receiver  of 
proper  dimensions  and  registered  by  an  incoherent  autodyne  detector.  An  analytical 
framework  is  developed  to  establish  the  fundamental  limits  on  the  resolution  of  such  a  system 
to  obtain  a  quantitative  understanding  of  the  dependence  of  the  accuracy  of  estimation  on 
target  separations,  amplitudes,  relative  velocities,  and  on  the  number  of  photons  available. 
The  value  of  a  priori  knowledge  to  obtain  a  desired  accuracy  for  estimation  of  a  parameter  of 
interest  is  computed  and  discussed  for  selected  cases.  The  analysis  is  carried  for  Poisson 
photon  statistics.  The  performance  of  various  estimators  is  evaluated  by  utilizing  the  Cramer- 
Rao  bound*®’**  (CRB)  for  both  shot  noise  and  background  limited  conditions.  The  results 
presented  have  applications  for  robust  laser  radars  and  contribute  to  a  better  understanding  of 
the  performance  and  potential  of  an  autodyne  detection  system. 

There  are  several  special  cases  corresponding  to  interesting  problems.  When  the 
relative  velocity  of  two  points  is  known  and  equal  to  zero  then  we  have  to  determine 
longitudinal  separation  of  two  fixed  points,  i.e.,  time  delay  between  echo  signals  registered  by 
the  detector.  This  problem  is  a  one  dimensional  temporal  analogue  of  the  estimation  of  a 
spatial  separation  of  two  point  sources  from  their  image  *2.13  Another  problem  corresponds 
to  the  case  when  separation  between  two  points  is  zero  and  when  Doppler  shift  gives 
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information  about  relative  velocity  of  the  target  component.  This  can  be  the  case  for  instance 
in  some  geometries  when  the  nozzle  and  reflection  plume  are  illuminated  by  a  laser  pulse. 
The  case  when  reflectivity  of  one  of  tne  point  targets  is  known  and  much  higher  than  another 
brings  autodyne  detection  into  analogy  with  heterodyne  detection. 

Analysis  presented  in  this  paper  will  remain  valid  if  the  assumptions  that  a  target 
consists  of  two  points  are  relaxed.  Indeed,  two  targets  illuminated  by  a  relatively  narrow  laser 
beam  can  have  finite  dimensions  in  the  directions  perpendicular  to  the  line  of  sight  (LOS)  of 
illumination  beam  (provided  that  they  do  not  rotate).  The  longitudinal  dimensions  of  the 
target,  however,  should  be  sufficiently  small  so  that  the  echo  signal  produced  by  the 
individual  target  component,  pulse  spreading  or  distortion  of  phase  modulation  are  negligible. 
Transverse  dimensions  of  a  two  component  target  will  affect  performance  of  an  autodyne 
detector  through  limitations  imposed  on  the  size  of  the  photon  bucket.  Indeed,  since  the 
photon-bucket  area  should  not  exceed  (X^/A-p^)R^  where  X-  wavelength,  R  -  range  and 
Af£-  transverse  area  of  illuminated  target  assuming  that  all  separations  are  much  smaller  than 
R)  the  transverse  target  dimensions  limit  the  power  collection  capability  of  individual  photon- 
buckets. 


One  of  the  important  problems  addressed  is  how  the  estimation  accuracy  of  a  subset 
of  parameters  depends  on  the  knowledge  (and  the  numerical  value)  of  the  remaining 
parameters.  Indeed,  it  is  important  to  recognize  that  the  longitudinal  resolution  of  two  point 
targets  (their  separation  and  relative  velocity)  depends  not  only  on  the  temporal  pulse 
distribution,  total  signal  and  background  noise  level  but  equally  important,  on  the  state  of  our 
knowledge  about  the  object  prior  to  measurement. 

Our  goal  is  to  determine  theoretical  limits  on  the  accuracy  of  estimating  the  values  of 
a  set  of  parameters  describing  the  echo  signals  from  a  two  point  source.  Analysis  of  the 
accuracy  with  which  one  can  measure  separations,  individual  positions  and  relative  velocity 
can  be  extended  to  other  parameters  or  their  combinations.  A  way  to  determine  the  number 


C49-701 


-46- 


RIVERSIDE  RESEARCH  INSTITUTE 


of  signal  photons  necessary  to  measure  echo  signal  parameters  with  desired  accuracy  is 
indicated.  We  assume  that  the  echo  signal  representing  backscattering  from  two  point  targets 
and  which  may  be  embedded  in  a  uniform  background  is  registered  by  an  incoherent  (square- 
law)  detector.  Calculations  are  performed  in  the  limiting  case  of  ideal  detectors  with 
negligibly  small  (compared  to  the  fine  structure  of  the  backscatter  signal)  response  time, 
detecting  the  entire  echo  signal.  The  analysis  is  carried  out  for  Poisson  photon  statistics.  The 
Cramer-Rao  bounds  (CRB)  for  various  parameters  describing  echo  signal  are  computed  and 
the  sensitivity  of  the  CRB  to  different  parameters  and  changes  in  the  number  of  unknown 
parameters  is  examined.  The  results  presented  here  provide  a  more  comprehensive  definition 
of  resolution  of  longitudinal  separation  and  relative  velocity  for  autodyne  signals  and  are 
important  in  such  applications  as:  astronomy,  light  communications  and  reconstruction  of 
images  from  the  recorded  spatio-temporal  evolution  of  speckle  patterns. 

3.2.1  Statistics  and  Parametrization  of  Echo  Signals 

Let  P(t,(A))  be  the  power  distribution  on  the  detector  as  a  function  of  time,  where  (A) 
=  (Aj,...,An)  denotes  a  set  of  deterministic  parameters  describing  the  echo  signal.  The 
normalization  of  P(t,(A))  is  such  that  the  probability  of  nk  photocounts  in  the  k-th  time  slot  on 
a  detector  with  integration  time  r  is  given  by: 

p(nkl(A)>  -  Jy  [Ek(A)J  k  •’  ,  (3.15) 

where:  k 

*k((-))  "  }  dt  P<*'  W>  *  Pr  •  Hk((A))  +  Nb  ,  (3.16) 

<T)k 

P  denotes  the  signal  and  p-uniform  background  power.  In  our  analysis  we  assume  that  the 
backscattered  signal  is  embedded  in  Poisson  background  noise  and  is  continuously  sampled 
through  the  entire  echo  signal.  We  assume  the  detector  response  time  much  smaller  th»n 
integration  time  of  the  detector  and  small  enough  so  we  can  write 
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*k((A))  -  r[P<tk, (A))  +  p]  -Nk((A))  +  Nb  (3.17) 

where  denotes  sampling  time  of  the  detector.  The  joint  probability  distribution  of  collected 
photocounts  is  given  by 


P(n I  (A)) 


Jc«l 


P(nkl (A)) 


(3.18) 


where  n  is  a  vector  representing  the  photocounts  registered  during  observation  time. 


Let  the  time  dependent  part  of  the  complex  amplitude  of  the  illuminating  pulse  be 
described  by  f(t): 


f  (t) 


Z(t)« 


itwQ7 (t) 


(3.19) 


where  Z(t),  ?(t)  describe  the  amplitude  and  fret,  .ency  modulation.  We  selected  Gaussian 
amplitude  and  linear  frequency  modulation: 


Z(t) 


2  2 
-  t  Ho 


i  (t)  -  (1  +  at) 


(3.20) 


Assuming  that  Doppler  shifts  are  due  to  the  linear  (and  nonrclativistic)  motion  of  two  target 
components,  the  power  intercepted  by  the  receiver  can  be  parameterized  as  follows; 


P (t, (A) )  -  NalZ2(t  -  r ^  +  N,2z2(t  '  T2) 

+  2-lN#1N#2  Z(t-r1)Z(t-r2)co»[2a>oa(r-r2)  (t-(rx+r2)/2) 

+  -  (rx+  r2)/2))  +  *] 

-  Na|pZ2(t-r-Ar/2)  ♦  (1-p) Z2 (t-r+Ar/2) 

+  2^p (1-p)  Z (t-r-dr/2) Z (t-r+Ar/2)  • 
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co«^Au(t-T)  ♦  2«o»Ar 
where  the  Doppler  shift  of  carrier  angular  frequency  is  given  by 

Ll>  -  «o(2Av/c)  -  2»[S~]  (3.22) 

with  Av  representing  the  relative  velocity  of  the  two  target  components  along  the  LOS;  Nsj  « 
NsP  (N$2  =  Ns(l-p))  represent  the  average  number  of  photons  backscattered  from  the  first 
(second)  target  registered  by  the  detector  when  another  target  is  absent;  r  j  and  t2  are  round 
trip  time  delays  and  ♦  is  the  phase  offset;  p  is  related  to  the  relative  backscattered  signal  from 
the  target  components,  r  is  round  trip  time  delay  r  =  (r  i+r^ri  of  the  transmitted  pulse  to 
the  midpoint  of  the  line  segment  connecting  two  points  along  the  LOS,  and  Ar  =  (r  j  -  rj)  is 
differential  time  delay  which  we  can  call  temporal  separation  of  two  point  targets.  In  deriving 
this  equation  we  neglected  several  terms  which  were  of  the  order  of  v-/c  (vj  -  velocities)  as 
compared  to  those  already  present.  Another  term  which  was  omitted  in  the  argument  of  the 
cos(  • )  function  was  of  the  form  aAw(t-r)2;  dropping  this  term  is  justified  provided  that 
ac«l  which  we  assume  is  satisfied.  One  may  add  that  although  the  definition  of  pulse 
duration  can  vary  depending  on  the  application,  its  values  typically  are  between  2 a  and  4 o. 


(t-r)  +♦ 


]} 


(3.21) 


The  set  of  parameters  (A)  describing  backscattered  signal  registered  by  the  detector 
consist  of 

A  ■  ( r.  At,  Ld,  p,  Pi  *  (A^,  Aj,  *  *  *  (3  *23) 

and  the  total  (average  number  of  signal  photons  registered  by  detector)  is  given  by 

■  Hgj1  +  2p'lp(l-p)  •*co*4,| 


- a  [Au  ♦  2u  aArJ2 

6  a 


(3.24) 


C49-701 


-49- 


RIVERSIDE  RESEARCH  INSTITUTE 


From  this  expression  we  can  see  that  Ns  is  equal  to  the  average  number  of  photons  only  when 

♦  =  0°  or  180°  and  approximate  equality  holds  when  separation  between  two  point  sources 
(along  the  LOS)  or  <r^[A«  +  2<u0aAr]^  is  sufficiently  large.  Out  of  six  parameters  describing 
the  pattern  of  echo  signal,  two  (r  and  N$)  define  the  pattern  translation  in  time  and  scaling, 
respectively,  and  ;l<e  remaining  four  define  the  pattern  shapes. 

In  some  calculations  under  shot  noise  limited  conditions,  one  has  to  proceed  with  a 
special  care  if  the  echo  pattern  has  zeros.  Such  zeros  of  pattern  are  present  (excluding  zeros 
for  t-.+«),  when  pattern  parameters  satisfy  the  relation: 

2 

(Aw+2a»o»Ar)  j^-  ln(-“]  J  -  (2k+l)ir  (3.25) 

where  k  is  an  arbitrary  integer  and  Ar  +  0. 

Various  cases  of  echo  signal  patterns  are  shown  in  Figs.  3.1,  3.2  and  3.3  for  the 
special  case  of  p  ■  1/2,  when  the  backscatter  signal  from  two  target  components  is  of  equal 
strength  and  in  the  absence  of  phase  modulation.  In  these  and  other  figures,  td,  tdd  and  vd 
denoted  time  delay  to  the  midpoint,  temporal  separation  of  two  point  targets,  and  relative 
angular  Doppler  shift  measured  in  units  of  a,  i.e.;  td  •  r/a,  tdd  ■  A r/a,  vd  *  A ua.  The  three 
special  cases  of  received  echo  signals  when  phase  ♦  »  0°  (solid  line)  t  ■  90°  (dashed  line) 
and  ♦  »  180°  (dotted  line)  are  shown  in  several  figures  and  this  convention  is  preserved 
throughout  this  work  unless  indicated  otherwise. 

In  Fig.  3.1  wc  displayed  the  set  of  patterns  obtained  when  relative  LOS  velocity  of 
two  points  is  zero,  i.e.,  relative  longitudinal  position  is  fixed  and  temporal  separation  varies. 
This  case  is  similar  to  the  problem  of  spatial  resolution^.  In  Fig.  3.1a  the  echo  pattern  for  ♦ 

*  180°  corresponds  to  zero  signal.  Echo  signal  for  ♦  -  0°  is  twice  the  signal  for  *  ■  90°.  In 
Fig.  3.1b  the  pattern  for  ♦  ■  180°  is  barely  discernible.  In  Fig.  3.1c  and  3.  Id  the  echo  signals 
corresponding  to  three  different  relative  phases  are  quite  distinct.  When  temporal  separation 
increases,  echo  signals  become  less  sensitive  to  the  value  of  phase,  and  in  Fig.  3. If 
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Fig.  3.1  Echo  signal  patterns  from  two  equal  strength,  stationary  backscatters,  with 
changing  values  of  separation.  Solid  line  phase  angle  *  *  90°;  dotted  line  • 
phase  angle  +  -  180°  for  this  and  all  subsequent  figures  (3.1  •  3.13). 
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Tint  (twit*  of  o ) 


Echo  signal  patterns  from  two  equal  strength,  moving  backscatters  with  changing 
values  of  relative  velocity  and  separation  value  of  zero. 
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Fig.  3.2 
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all  three  patterns  are  almost  identical,  as  the  overlap  between  signals  from  two  components  of 
the  target  becomes  negligible.  Fig.  3.2  shows  the  echo  patterns  for  maximum  overlap  of 
backscattered  fields  when  temporal  separation  of  targets  is  zero  and  for  different  relative 
velocities.  Note  the  similarity  of  patterns  corresponding  to  ♦  a  180°  in  Fig.  3.2b  and  3.1c 
with  the  pattern  corresponding  to  the  same  phase  in  Fig.  3.1b  and  3.1c.  In  Fig.  3. 2d  we 
observe  asymmetry  for  the  case  corresponding  to  *  =  90°.  As  relative  velocity  increases, 
echo  signals  are  becoming  more  oscillatory,  the  width  of  individual  peaks  narrows  (— 1/vd) 
and  the  envelopes  of  oscillatory  pattern  are  provided  by  patterns  for  vd  =  0,  ♦  =  0°  and  *  * 
180°  (for  arbitrary  p). 

In  Fig.  3.3  we  plotted  echo  signals  for  selected  cases  of  nonzero  temporal  separations 
and  relative  angular  Doppler  shifts.  While  some  patterns  are  quite  different  from  those  in  Fig. 
3.1  and  3.2,  some  are  similar.  For  instance,  comparing  Fig.  3.1e  and  Fig.  3.3f  we  expect  it 
may  be  difficult  to  make  proper  pattern  classification  in  the  presence  of  noise.  As  we  could 
expect,  only  patterns  corresponding  to  ♦  =  180°  have  zeros.  It  is  important  to  note  that 
patterns  with  p  »  1/2  have  a  structure  that  is  most  pronounced,  for  p  approaching  0  (or  1)  the 
echo  signal  will  be  represented  by  a  Gaussian  shape  with  weak  modulations. 

3.2.2  Performance  Measures 


The  performance  of  the  estimators  of  pattern  parameters  is  evaluated  by  utilizing  the 
Cramer-Rao  bounds^.  These  bounds  are  simple  lower  bounds  on  variances  of  unbiased 
estimators.  In  this  report  we  study  the  sensitivity  of  the  CRB’s  to  the  temporal  separation  of 
point  targets,  their  amplitudes,  relative  phase  and  relative  velocity  (Doppler  shift). 

In  general,  if  parameters  (A)  are  arranged  in  such  a  way  that  the  first  m  components 
of  vector  A  denote  unknown  deterministic  parameters,  then  the  Cramer-Rao  bounds  on  the 

A 

variance  of  unbiased  estimators  Aj(i  *  l,...,m)  are  computed  from  the  m  x  m  Fisher 
information  matrix  J((A))  with  elements  J((A))jj  given  by: 
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. .  .  /llnp(nl(A))  31np(nl(A))\ 

-  <  N - ^  3^  /n 


(3.26) 


OAj 

where  <  *>n  stands  for  averaging  over  all  realizations  of  n  governed  by  the  distribution 
p(n  I  (A)). 


The  CRB’s  for  variance  of  estimators  A  are  given  by  the  diagonal  elements  of  the 
inverted  information  matrix  denoted  as  (r^((A)))jj  so  that: 


r-1 


var  i  ( J  A ( (A) ) ) -  CRB (A^ (A) ) 


(3.27) 


Set  (A)  of  deterministic  parameters  can  contain  r,  A r,  Ao,  p,  N$  and  ^  (or  related 
parameters).  For  the  probal  .lity  distribution  functions  (PDF)  p(n  I  A)  which  are  the  products 
of  PDF’s  such  as  those  in  **).  (3.18),  the  Fisher  information  matrix  J((A))jj  is  given  by: 

E:  (alnP(nC>-31PP,°3xf->  >n  •  (3  28> 
1c*  1.  •  J 

When  signal  statistics  arc  governed  by  Poisson  distribution,  Eq.  (3.28)  becomes: 

L 


J((A)) 


ij 


1  3<nv>  3<nw> 


fei  <n&  3Ai 


3  A. 


j 


(3.29) 


where  <n^>  denotes  the  expected  values  of  registered  photocounts. 


In  the  limit  of  small  integration  times  covering  an  entire  echo  signal,  Eq.  (3.29)  can  be 
written  in  the  following  form: 


J((A,)ij  -  | 


dt 


3 

3 

3*7  r(t'  <*>> 

nr  *»•«*» 

i 

j 

P(t,  (A) )  +p 


(3.30) 
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This  expression  remains  valid  when  the  observation  time  is  finite,  provided  that  infinite 
time  domain  is  replaced  by  integration  over  a  finite  one. 


The  selection  of  parameters  describing  the  pattern  of  interest  depends  on  the 
application.  In  general,  the  Fisher  information  matrix  Jy  for  parameters  Bj  =  Bj  ((A))  (i  =  1, 
....  m)  is  related  to  the  Fisher  information  matrix  Jy  for  parameters  Aj  (i  *  l,...,m)  as  follows: 


J((A)) 


ij 


k,l 


9B  ((A))  _ 

5x7“ 


kl 


( (A) ) 
9a" 


(3.31) 


Two  important  special  cases  are  obtained  when  the  power  density  of  the  uniform 
background  is  negligible  compared  to  the  power  of  the  echo  pattern  and  when  the  background 
power  is  much  larger  than  the  signal  power.  We  will  refer  to  these  two  limiting  cases  as 
signal  shot  noise  and  background  shot  noise  limited,  respectively.  In  the  second  case,  Eq. 
(3.30)  for  the  Fisher  information  matrix  may  be  approximated  by: 


J((A))i;J« 


\  \ 


dt 


3A. 


P(t,  ((A)) 


3A 


P(t,  (A)) 


j 


(3.32) 


i-  )  "n 


where  we  assume  that  the  set  of  estimated  parameters  does  not  include  p. 


3.2  J  Expressions  Used  In  Numerical  Calculations 


When  computing  the  Fisher  information  matrix,  for  the  most  general  case  of  six 
unknown  parameters,  twenty-one  integrals  have  to  be  computed  since  matrix  J((A))y  is 
symmetric.  Each  integral  may  depend  on  up  to  seven  parameters  (including  background 
density  p).  Diagonal  matrix  elements  of  the  inverse  of  this  6x6  matrix  provide  the  bounds 
on  parameter  estimation.  Note  that  for  small  temporal  or  small  velocities  separations  of  point 
targets,  the  Fisher  information  matrix  can  become  ill-conditioned. 

In  order  to  achieve  numerical  stability,  it  may  bt  helpful  to  find  the  relationships 
between  various  matrix  elements.  To  reduce  the  amount  of  numerical  computations,  and  to 
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clarify  the  dependence  on  various  parameters  defining  the  pattern  of  interest,  it  is  useful  to 
identify  those  parameters  which  can  be  removed  by  appropriate  transformations  of  matrix 
elements.  It  is  also  useful  to  determine  the  minimum  number  of  significant  parameters. 
Examining  Eq.  (3.30),  we  note  that  by  multiplying  elements  of  the  Fisher  information  matrix 
by  some  power  of  Nso,  we  can  make  these  elements  dependent  on  Ns0  only  through  a  factor 
of  p/N^.  This  scaling  property  allows  the  following  factorization  of  the  bounds: 

N  CRB  (A.)  *  N  £%  [  r,  At,  Cm,  p,  N  ,  *,  p 

9  i  »  Aj  (  S 

Ar,Aw,  p,  Ns'*/P/Na]  (3.33») 


for  Aj»*  N$,  and 


Hr- <*»(",)  -  Hr  vH 

-  fM  [r,Ar,Au,p,  1,  *,  p/N^j  (3.33b) 

If  the  observation  time  covers  the  entire  echo  pattern,  in  the  limit  of  small  integration 
times,  we  observe  that  all  matrix  elements  are  r  independent.  This  is  due  to  the  fact  that  the 
dependence  on  r  in  Eq.  (3.30)  can  be  removed  by  a  simple  change  of  variables  (consisting  of 
translation  by  r).  Without  any  loss  of  generality,  we  may  set  r  equal  to  0  in  Eqs.  (3.33a)  and 
(3.33b). 


These  considerations  lead  to  the  reduction  of  the  number  of  nontrivial  parameters 
from  seven  to  five,  i.e.,  Ar,  Au,  p,  ♦,  p/Ny  In  the  shot  noise  limited  case,  the  number  of 
parameters  is  equal  to  four.  This  is  also  the  number  of  parameters  in  background  limited  case 
if  the  CRB’s  are  scaled  by  an  additional  factor  of  Ng/p  (or  Ng/Ng  where  Ng  =  pa  if  we  adopt 
our  convention  for  time  scale)  to  obtain  results  that  are  independent  of  Ng/p. 

Since  the  Fisher  information  matrix  can  be  ill-conditioned,  it  is  useful  to  establish 
additional  relations  in  special  cases  in  order  to  detect  possible  numerical  instabilities  which 
may  occur.  An  important  special  case  corresponds  to  the  signal  shot  noise  limited  conditions. 
If  we  set  p  =  0,  the  following  matrix  elements  can  be  easily  evaluated  (see  Eq.  (3.30): 
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< 


(3.34) 


In  addition,  an  interesting  special  case  under  signal  shot  noise  limited  conditions  is 
obtained  if  we  assume  the  value  of  *  =  0  or  t  (*  maybe  known  or  unknown)  and  the 
remaining  part  of  the  argument  of  cos(  • )  function  in  Eq.  (3.21)  identically  vanish,  which 
happens  when  Aw  +  awQAr  =  0.  In  this  case  most  of  the  matrix  elements  of  the  Fisher 
information  matrix  can  easily  be  evaluated.  Indeed,  let  us  assume  for  simplicity  that  Aw  =  0, 
then,  all  nontrivial  matrix  elements  Jjj;  Aj,  Aj  *  ♦,  Aw  can  be  explicitly  evaluated. 


P(t;(A))  -  ((p  Z(t  ♦  r  -  |r)+  H(l-p)Z(t  ♦  r  ♦  ^~)  ]  } 

-  (L (t ; (A) ) ) 2  ;  e  ■  cos*  (3.35) 

Therefore,  the  information  matrix  J((A))jj  (Aj,  Aj  *  Aw,*)  can  be  written  in  the  following 
form: 


J((A)) 


Vj 


CD 

|  dt  ||-L(t;(A)) 

£r-  L(t;<A>) 

l5Xj  '  1 

(3.36) 


-  4  -§ - —  [*t  Mt;(A))  L(t;  (A) ) 

3Ai  3A^  j 


(A)  -  (A) 


Since  L  is  the  sum  of  two  Gaussians  and  the  product  of  two  L-functions  is  again  a 
combination  of  Gaussians,  the  integral  in  Eq.  (3.36)  can  be  easily  evaluated.  After  taking 
appropriate  derivatives,  we  obtain: 


JN  N  -  (1/H#>  <1  ♦  2r  T(Ar,p)),  JN  -  0 
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where 


T (Ar, p)  *  ’lp(l-p)  sxp  -(Ar)  /8a2 


(3.37b) 


Another  special  case  in  the  shot  noise  limited  regime  is  obtained  when  Ar  =  0  and  p  -  1/2. 
Again,  one  can  write  P(t,(A»  as  a  square  of  the  simple  expression 
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2N  ‘1/4 
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(3.38) 


and  the  information  matrix  can  be  written  in  a  form  analogous  to  Eq.  (3.^6)  wit-,  tunction  l 
replaced  by  L\  After  evaluating  the  integral  and  its  derivatives,  we  obtain: 
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JtAw"  ”  Ns  s<Aw)ain*,  "  N,[1  "  S(Aq)co#*]  , 

m  S  (&»>)  ain#,  3 *  "Au  (3.39a) 

where 

S(Au)  ■  axp(- (aA«)2/2)  (3.39b) 

Similar  formulas  can  be  obtained  for  the  background  limited  case  when  the  information 
matrix  J((A))jj  can  be  written  in  the  following  form: 


J (  (A)  )  -  (1/p) - 

‘  Vj  3A. 


_3 

3A 


--  Jp(t,  (A)  )  P  (t ,  (A))dt 

j  ~oo  A 


(3.40) 


Again,  the  integral  can  be  evaluated  and  derivatives  performed  leading  to  closed  form 
expressions  for  all  matrix  elements.  These  important  formulas  are  more  complicated  and  will 
be  presented  and  discussed  elsewhere. 


In  numerical  calculations  of  the  Fisher  information  matrix,  an  integration  time 
interval  of  15  a  was  used.  The  integration  grid  consisted  of  1000  points.  The  loss  of  accuracy 
due  to  a  finite  integration  region  and  coarseness  of  the  integration  grid  did  not  exceed  a  few 
percent  for  the  considered  values  of  the  parameters. 

In  some  trials  for  very  small  temporal  separations  and  velocities  (tddiO.l,  vdSO.l) 
the  Fisher  information  matrix  became  numerically  unstable  and  quadruple  precision  was 
necessary  to  obtain  meaningful  results. 


3.2.4  Discussion  of  Results 


The  sensitivity  of  estimation  for  any  given  parameter  depends  on  the  actual  v  ?'.«<  *  o 
that  parameter,  and  on  parameters  describing  a  two  point  target.  In  additio  v  ;e  curacy  jf 
estimation  depends  on  whether  the  other  parameters  are  known.  Or...  rc.  F  iiiei  information 
matrix  is  computed  for  the  general  case,  the  bound'  can  v  c  ootainrc  ior  special  cases  by 
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inverting  an  appropriate  submatrix  when  some  of  the  parameters  are  known.  Eq.  (3.31) 
allows  us  to  derive  the  transformed  Fisher  information  matrix  for  any  function  of  our  selected 
canonical  parameter  set  and  thus  obtain  the  desired  bounds. 

Two  limiting  cases  discussed  here  correspond  to  the  signal  shot  noise  and 
background  limited  cases.  In  these  two  limits  the  signal  dependence  is  easily  scaled  out  and 
results  depend  on  the  values  of  only  four  parameters:  temporal  separation  of  two  targets  dr, 
fractional  strength  of  the  first  target  P,  and  the  relative  phase  angle  ♦  and  angular  Doppler 
frequency  shift  bu. 

The  following  figures  describe  the  accuracy  of  estimation  of  the  echo  for  various 
pattern  parameters  as  a  function  of  temporal  separation  between  two  point  targets  and  relative 
Doppler  shift.  Six  consecutive  figures,  3.4  through  3.9,  are  arranged  in  similar  fashion:  the 
left  panels  refer  to  cases  when  all  parameters  are  known  except  the  one  being  estimated;  the 
right  panels  refer  to  cases  when  all  parameters  are  unknown.  All  panels  correspond  to  the  true 
value  of  p  =  1/2  and  each  panel  contains  three  graphs  corresponding  to  different  values  of 
relative  phase;  *  *  0°  (solid  line)  +  *  90°  (dashed  line)  and  *  «  1 80°  (dotted  line)  used 
before.  Alternate  sets  of  figures  show  signal  shot  noise  limited  cases  and  background  shot 
noise  limited  cases,  respectively.  Note  that  dependence  on  signal  value  was  scaled  out.  In 
these  figures,  Ng  denotes  the  average  number  of  background  photons  registered  in  time  a. 

We  first  consider  the  estimator  for  the  temporal  separation  of  two  point  targets.  A 
comparison  of  three  cases  in  Fig.  3.4a  shows  behavior  which  may  be  surprising:  for  small 
separations  it  is  advantageous  to  have  point  targets  in  phase  opposition  rather  than  in  phase  or 
phase  quadrature.  This  is  despite  the  much  smaller  number  of  photons  available  in  the  pattern 
(see  Fig.  3.1a  -  3.1b).  The  apparent  contradiction  is  easily  explained  since  this  is  a  case  in 
which  all  other  parameters  and,  in  particular,  the  expected  number  of  photons  are  known.  It  is 
this  small  number  of  registered  photons  that  indicates  that  the  temporal  distance  between  point 
targets  must  be  small.  Indeed,  in  Fig.  3.4b,  corresponding  to  the  case  in  which  all  parameters 
are  unknown,  the  performance  of  the  estimator  of  temporal  separation  is  reversed;  it  is  best 
when  targets  are  in  phase,  and  worst  when  they  are  in  phase  opposition. 
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Fig.  3.4  Dependence  of  scaled  CRB  (tdd)  on  temporal  separation  of  two  point  targets  of 
equal  strength;  selected  cases  for  different  values  of  relative  velocity.  Signal  shot 
noise  limited  conditions. 
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Fig.  3.5 
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Dependence  of  scaled  CRB  (tdd)  on  temporal  separation  cf  two  point  targ 
equal  strength;  selected  ca^s  for  different  values  of  relative  vel 
Background  shot  noise  limited  conditio  s. 
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Dependence  of  scaled  CRB  (v<l)  on  temporal  separation  of  two  point  targets  of 
equal  strength;  selected  cases  for  different  va)a;.s  of  relative  velocity.  Signal  shot 
noise  limited  conditions. 
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Fig.  3.7  Dependence  of  scaled  CRB  (vd)  on  temporal  separation  of  two  point  targets  of 
equal  strength;  selected  cases  for  different  values  of  relative  velocity. 
Background  shot  noise  limited  conditions. 
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Fig.  3.8  Dependence  of  scaled  CRB  (tel)  on  temporal  separation  of  two  point  targets  of 
equal  strength;  selected  cases  for  different  values  of  relative  velocity.  Signal  shot 
noise  limited  conditions. 
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Fig.  3.9  Dependence  of  scaled  CRB  (td)  on  temporal  separation  of  two  point  targets  of 
equal  strength;  selected  cases  for  different  values  of  relative  velocity. 
Background  shot  noise  limited  conditions. 
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Fig.  3.5a  (background  limited  case),  shows  that  for  small  separations,  performance  is 
best  for  the  case  when  targets  are  in  phase,  and  slightly  worse  when  the  targets  are  in  phase 
opposition  or  quadrature.  This  is  what  one  would  expect;  in  the  presence  of  background 
generated  shot  noise,  absence  of  the  signal  is  not  a  very  good  indicator  of  temporal  separation. 
Figs.  3.4c,  3.4e,  and  3.4g  describe  the  performance  of  estimator  of  temporal  separation  for 
nonzero  relative  target  velocity.  In  contrast  to  the  static  case,  selection  of  ♦  »  180°  no  longer 
gives  the  best  performance  of  the  estimator.  The  presence  of  relative  motion  can  lead  to 
improved  performance  for  some  values  of  temporal  separations  and  phases,  and  degraded 
performance  for  others.  It  is  important  to  note  that  for  angular  Doppler  shifts  (vd)  of  the  order 
of  four  and  larger,  the  performance  of  estimators  of  temporal  separations  becomes  phase  and 
velocity  independent  (Fig.  3.4g).  Since  for  large  temporal  separation  overlap  of  backscattered 
fields  from  two  target  components  is  negligible,  independence  of  the  estimator  of  separation 
on  phase,  velocity  and  temporal  separation  is  expected  for  large  tdd.  This  behavior  is  indeed 
observed  in  all  panels  of  Figs.  3.4  and  3.5. 

When  we  examine  Fig.  3.4b,  3.4d,  3.4f  and  3.4h,  where  all  parameters  are  unknown 
and  compare  them  with  the  corresponding  cases  when  all  parameters  but  separation  are 
known,  we  will  notice  that  strong  sensitivity  to  phase  decreases  with  increasing  relative 
velocity  and  the  very  poor  estimator  performance  (for  small  relative  velocities)  improves  as 
well.  Indeed,  when  vd  =  4.0,  Fig.  3.4g  and  3.4h  are  practically  identical,  showing  little  value 
of  a  priori  knowledge  in  this  case.  This  is  in  contrast  with  the  static  case,  when  two  sources 
arc  in  quadrature  (#  =  90°)  and  temporal  separation  is  1  a;  a  comparison  of  Fig.  3.4a  and  3.4b 
shows  that  we  need  approximately  twenty-five  times  more  photons  to  compensate  for  the  lack 
of  a  priori  knowledge  of  the  other  parameter  values. 

Computation  of  the  minimum  average  number  of  registered  photons  necessary  to 
perform  a  measurement  of  temporal  separation  with  given  accuracy  is  straightforward  when 
all  other  parameters  are  known.  For  instance,  let  p  =  0.5,  vd  =  4.0,  *90°,  (so  that  CRB  is 

described  by  Fig.  3.4g  and  assume  that  we  want  to  measure  temporal  sepo.  .ions  of  tdd  *  0.5 
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with  accuracy  0.1  (l<r  error).  From  Fig.  3.4g  we  can  determine  that  N$CRB(tdd)  is 

A 

approximately  10  for  temporal  separation  O.S.  Since  the  CRB  (tdd)  has  to  be  smaller  than 
(O.l)^  m  qi  for  tj,e  required  accuracy,  the  minimum  required  number  Ns  of  photons  is  1000. 
When  more  parameters  are  unknown,  such  calculations  may  not  be  sufficient  and  may  require 

A 

for  instance  finding  the  bound  on  CRB  (tdd)  when  other  parameters  vary  within  the  required 
range.  (See  some  additional  comments  in  the  discussion  of  midpoint  time  delay  estimation). 

Note  that  for  very  small  separation,  the  estimation  bounds  in  some  cases  seem  to 
diverge.  In  Fig.  3.4a,  we  can  see  that  for  small  separations  and  sources  in  phase  opposition, 
the  bound  has  a  finite  limit,  while  for  sources  in  phase,  we  observe  divergent  behavior.  This 
also  seems  to  be  true  for  sources  in  quadrature.  One  can  analytically  verify  this  behavior.  For 
the  case  when  all  parameters  are  known,  divergent  behavior  is  often  observed  only  for  special 
parameter  values,  e.g.,  p  =  1/2.  When  all  parameters  are  unknown,  however,  divergent 
behavior  for  small  separations  is  prevalent.  This  divergence  is  related  to  the  nonexistence  of 
unbiased  parameter  estimators  for  small  temporal  separations.  This  behavior  for  small 
separations  will  be  discussed  elsewhere.  The  present  remarks  apply  as  well  to  the  behavior 
observed  for  the  other  parameter  estimators. 

With  the  exceptions  mentioned  already,  Fig.  3.4  and  Fig.  3.5  are  remarkably  similar. 
In  the  background  limited  case  temporal  separation  estimators  seem  to  be  more  phase 
sensitive  and  their  variance  may  vary  within  a  larger  range. 

The  next  estimator  considered  is  for  angular  Doppler  frequency  shift  shown  in  Fig. 
3.6  and  3.7  for  signal  shot  noise  and  background  limited  conditions,  respectively.  The  first 
observation  is  that  Fig.  3.6a  as  well  as  Fig.  3.6b,  3.7a  and  3.7b  display  only  one  bound 
corresponding  to  ♦  *  90°.  This  is  because  finite  variance  unbiased  estimators  of  vd  do  not 
exist  for  ♦  =  0°  and  180°  and  for  relative  velocities  approaching  zero.  It  is  not  difficult  to 
understand  why  this  is  the  case.  Angular  Doppler  frequency  shift  enters  echo  signal  in  Eq. 
3.38  through  a  cos(.)  factor.  For  phase  angle  t  ®  0°  and  180°  and  for  very  small  velocities, 
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expansion  of  cos(.)  leads  to  a  term  proportional  to  (A w)2.  For  4*  =  90°  in  similar  conditions, 
the  first  term  in  expansion  is  proportional  to  Aw.  This  indicates  that  an  unbiased  estimator  of 
Au  may  exist  for  *  =90°  but  not  for  ♦  =  0°  or  190°  (and  an  unbiased  estimator  of  (Aw)2  may, 
however,  exist  for  these  two  phase  angles).  Already  for  small  relative  velocities,  however, 
corresponding  to  vd  =  0.5,  we  have  a  full  set  of  bounds  for  three  phase  angles.  In  fact,  in  the 
case  of  known  values  of  pattern  parameters,  the  dependence  of  the  Doppler  shift  estimator  on 
phase  angle  is  rather  weak  both  for  signal  shot  noise  (Fig.  3.6c,  3.6e,  3.6g)  and  background 
(Fig.  3.7c,  3.7e,  3.7g)  limited  performance.  For  relative  velocities  of  four  (Fig.  3.6g  and  3.7g) 
corresponding  to  vd  of  the  order  and  larger,  the  performance  of  the  estimator  of  relative 
Doppler  shift  becomes  phase  and  velocity  independent  in  full  analogy  with  the  performance 
of  the  estimator  of  temporal  separation.  The  dependence  on  phase  angle  is  much  stronger 
when  all  parameters  are  unknown  (Fig.  3.6d,  Fig.  3.6f,  Fig.  3.7d,  Fig.  and  3.7Q  but  this 
sensitivity  rapidly  diminishes  with  increasing  velocity,  and  for  vd  =  4  (Fig.  3.6h  and  3.7h)  the 
performance  of  an  estimator  of  relative  Doppler  shift  becomes  essentially  the  same  case  when 
knowledge  about  the  values  of  other  parameters  is  available.  As  we  could  expect,  the 
performance  of  estimator  of  relative  Doppler  shift  rapidly  deteriorates  with  the  increase  of 
temporal  separation  between  point  targets.  The  flatness  of  the  performance  curve  in  Fig.  3.6b 
and  3.7b  is  due  to  the  choice  of  performance  scale,  but  it  is  interesting  to  note  that  optimal 
performance  for  small  velocities  is  often  obtained  for  rather  large  temporal  separation  of  point 
targets. 


The  next  estimator  considered  is  for  midpoint  time  delay.  In  the  static  case  when  all 
parameters  are  known,  Fig.  3.8a  and  3.9a  show  the  expected  behavior  for  small  separations;  it 
is  more  advantageous  to  have  target  components  in  phase  rather  than  in  quadrature  or  phase 
opposition.  This  is  consistent  with  the  larger  number  of  photons  available  in  such  patterns. 
Different  behavior  is  observed  (for  small  separations)  when  all  parameters  are  unknown:  the 
performance  of  the  estimator  is  best  when  targets  are  in  phase  opposition  and  worst  when  they 
are  in  phase.  This  result  can  be  explained  by  correlation  of  this  estimate  with  an  estimator  of 
P  (see  discussion  in  reference  13).  As  expected,  for  large  separation,  the  performance  becomes 
phase  and  relative  velocity  independent  and  approaches  a  common  limit  in  all  :ases  whether 
parameters  are  known  or  unknown.  Rather  unexpected  and  exotic  behavior  is  observed  in 
cases  when  all  parameters  are  unknown  and  relative  velocity  is  nonzero.  In  both  signal  shot 
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noise  and  background  limited  cases  for  ♦  =180°,  the  performance  bound  shows  a  sharp  peak 
(Fig.  3.8d,  3.8f,  3.9d  and  3.9f).  This  behavior  is  quite  common  for  small  values  of  P  and  is 
due  to  various  exponential  and  trigonometric  factors  in  elements  of  the  Fisher  information 
matrix,  which  may  lead  to  rapid  change  of  behavior.  With  increase  of  velocity,  the 
dependence  on  phase  is  decreasing,  but  in  contrast  to  two  previous  estimators,  no  common 
limiting  performance  curve  is  attained  for  the  cases  when  pattern  parameters  are  known  or 
unknown.  In  fact,  one  can  show  that  for  small  temporal  separations  when  all  parameters  are 
known,  CRB  (tel)  decreases  as  l/(vd)^  (see  Eq.  (3.39a)  for  large  enough  velocities.  This 
supports  our  intuitive  feeling  that  it  is  easier  to  measure  position  displacement  of  a  more 
structured  pattern.  On  the  other  hand,  this  example  shows  that  some  care  is  necessary  in 
interpreting  CRB  predictions.  Indeed,  since  CRB  is  decreasing  like  l/(vd)^,  in  order  to 
perform  a  measurement  of  pattern  position  with  given  accuracy,  we  will  require  less  than  one 
photon  if  relative  Doppler  shift  is  sufficiendy  large.  This  is  obviously  not  a  sensible  answer. 
In  order  to  properly  use  the  prediction  given  by  CRB,  we  should  remember  that  it  is  a  local 
bound  useful  only  for  small  parameter  changes  (compared  to  fine  structure).  To  use  CRB  for 
prediction  of  pattern  shifts  in  the  present  application,  we  should  be  interested  in  shifts  with 
errors  smaller  than  the  distance  between  interference  fringes,  since  this  distance  is  of  the 
order  of  1/vd,  we  can  easily  establish  that  sufficiently  accurate  measurement  requires  a  few 
dozens  of  photons.  If  we  examine  behavior  of  the  midpoint  time  delay  estimator  in  the  case 
when  all  parameters  are  unknown,  behavior  of  CRB  for  small  displacement  seems  to  indicate 
that  it  is  also  decreasing  with  increase  of  relative  Doppler  shift.  However,  when  we  examined 
the  behavior  of  the  bound  for  vd  =  5  and  vd  ■  6  we  observed  no  further  significant  decrease 
from  the  value  observed  for  vd  *  4. 

Establishment  of  asymptotic  behavior  of  CRB’s  for  large  vd  is  certainly  of  interest . 
Indeed,  if  observed  stagnation  of  CRB  is  real,  arguments  similar  to  those  used  in  the  case 
when  all  parameters  are  known  would  indicate  that  the  number  of  photons  required  for  fine 
pattern  shift  estimation  increases  as  (vd). 


C49-701 


-71  - 


RIVERSIDE  RESEARCH  INSTITUTE 


Figures  3.10  and  3.1 1  are  arranged  in  different  fashion  than  the  previous  six  figures. 
All  panels  refer  to  the  cases  when  signal  is  unknown  and  p  =  0.01.  The  top  four  panels 
describe  the  performance  of  the  estimators  of  angular  relative  Doppler  shift  for  selected  values 
of  Doppler  shift  as  a  function  of  temporal  separation  between  point  targets,  while  the  bottom 
four  describe  the  performance  of  estimator  of  temporal  separation  for  similar  cases.  This 
important  case,  when  backscatter  signal  from  one  target  component  has  much  higher  strength 
than  return  from  the  other  component,  exhibits  highly  structured  behavior  of  performance 
curves.  Comparison  of  Fig.  3.10  for  signal  shot  noise  limited  case  and  Fig.  3.1 1  for 
background  limited  performance  shows  very  similar  behavior.  As  relative  velocity  increases 
in  a  similar  way  as  for  p  =  1/2,  the  performance  of  estimators  becomes  phase  and  relative 
velocity  independent  Similar  plots  for  the  case  when  all  parameters  are  known  (not  shown) 
show  less  dependence  on  phase  angle  than  for  p  =  0.5,  and  for  high  relative  velocity  they 
approach  the  same  phase  and  velocity  limit.  One  can  also  observe  a  correlation  between  some 
of  the  kinks  and  brakes  in  performance  curves  of  the  estimator  of  Doppler  and  separation.  A 
more  detailed  discussion  of  this  behavior  in  the  background  limited  case  will  be  presented  in  a 
forthcoming  paper,  when  analytic  results  will  be  available  for  the  Fisher  information  matrix. 

The  next  two  figures  illustrate  performance  of  estimators  of  relative  Doppler  shift 
and  separation  as  a  function  of  shot  noise  limited  performance.  Fig.  3.12  is  for  p  »  1/2,  and 
panels  on  the  left  are  for  signal  parameters  known  while  panels  on  the  right  are  for  unknown 
signal  parameters.  This  figure  provides  an  excellent  illustration  of  how  with  increase  of 
relative  velocity  the  performance  of  estimators  becomes  velocity  and  phase  independent. 
Similar  illustration  is  provided  in  Fig.  3.13  where  all  panels  arc  for  unknown  signal  and 
panels  on  the  left  are  for  p  *  0.1  and  on  the  right  for  p  *>  0.01.  Again,  one  can  see  a 
correlation  between  various  abrupt  changes  in  the  performance  of  various  estimators.  Another 
interesting  type  of  plots  of  estimator  performance  is  to  present  their  behavior  as  a  function  of 
phase  angle.  These  figures  are  not  shown  here  but  they  demonstrate  that  best  (or  worst) 
performance  corresponds  in  general  (for  vd  *  0)  to  phase  angles  different  than  three  standard 
phases  selected  in  our  figures. 
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Fig.  3. 10  Dependence  of  scaled  CRB  (vcl)  and  CRB  (tdd)  on  temporal  separation  of  two 
point  targets  of  differing  strength;  selected  cases  for  different  values  of  relative 
velocity.  Signal  shot  noise  limited  conditions. 
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Fig.  3.1 1  Dependence  of  scaled  CRB  (vd)  and  CRB  (t3d)  on  temporal  separation  of  two 
targets  of  differing  strength;  selected  cases  for  different  values  of  relative 
velocity.  Background  shot  noise  limited  conditions. 
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Fig.  3.12  Dependence  of  scaled  CRB  (vd)  and  CRB  (tdd)  on  relative  angular  Doppler  shift; 


of  two  point  targets  of  equal  strength;  selected  cases  for  Different  values  of 
temporal  separation.  Signal  shot  noise  limited  conditions. 
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at  tv.v  point  targets  of  different  strength;  selected  cases  for  different  values  of 


temporal  separation.  Signal  shot  noise  limited  conditions. 
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The  value  of  a  priori  knowledge  in  parameter  estimation  is  illustrated  in  Fig.  3.14. 
This  figure  shows  the  sensitivity  of  midpoint  time  delay  and  relative  angular  Doppler  shift 
estimation  for  varying  numbers  of  unknown  parameters  for  signal  shot  noise  limited 
performance.  Panels  on  the  left  are  for  phase  ¥  =  0°  and  on  the  right  for  t  =  180°.  Top 
panels  and  for  p  =  1/2  and  bottom  panels  are  for  p  =  0.01. 

When  we  examine  Fig.  3.14a  we  observe  that  knowledge  of  P  (cases  b)  or  ♦  (case  c) 
brings  substantial  improvement  in  estimator  performance  in  some  ranges  of  temporal 
separation.  When  both  parameters  arc  known  (case  d)  performance  is  essentially  the  same  as 
in  case  of  knowledge  of  all  parameters.  Fig.  3.4b  shows  a  more  dramatic  change.  Knowledge 
of  p  (case  a)  or  (case  b)  leads  to  the  performance  close  to  one  obtained  with  knowledge  of 
all  nuisance  parameters  (case  d).  In  fig.  3.14c  and  3.14d  we  observe  more  gradual 
improvement  with  increase  of  a  priori  knowledge.  This  improvement  can  be  expressed  in 
terms  of  smaller  number  of  photons  necessary  to  perform  the  measurement  with  required 
accuracy.  Such  considerations  provide  good  illustration  of  the  observation  that  achievable 
resolution  depends  not  only  on  signal  value  but  also  on  the  amount  of  a  priori  knowledge 
available  about  the  echo  pattern. 

3.2.5  Summary  and  Conclusions 

The  fundamental  limits  were  established  for  the  accuracy  of  the  estimation  of  range 
and  velocity  resolution  for  autodyne  detection  of  a  two  point  target  with  pulsed  laser  radar 
consisting  of  a  transmitter  collocated  with  a  photon  bucket  receiver.  The  Cramer-Rao  bounds 
(CRB’s)  were  computed  for  both  the  signal  shot  noise  and  background  limited  cases.  From 
these  bounds,  the  number  of  signal  photons  necessary  to  estimate  a  parameter  with  a  desired 
accuracy  may  be  determined.  The  study  of  longitudinal  position  estimation  can  be  of  interest 
in  tracking  applications. 

The  sensitivity  of  the  bounds  to  different  parameter  values  and  changes  in  the 
number  of  unknown  parameters  was  established.  The  value  of  a  priori  knowledge  in 
parameter  estimation  was  demonstrated  and  related  to  the  number  of  unknown 
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Fig.  3.14  Dependence  of  scaled  CRB  (td)  and  CRB  (v'd)  on  the  number  of  unknown 
parameters.  Panels  A  and  C:  phase  angle  *  ■  0°;  panels  B  and  D:  phase  angle  ♦ 
«  180°.  For  panels  A  and  B,  curve  designations  are:  a-  all  parameters  unknown; 
b  •  all  parameters  except  p  unknown;  c  -  all  parameters  except  t  unknown;  d  -  all 
parameters  except  p  and  9  unknown  (in  this  case  identical  curves  are  obtained  if 
all  parameters  are  known).  For  panels  C  and  D,  curve  designations  are:  a,  b,  c 
and  d  •  as  in  panels  A  and  B;  e  •  all  parameters  except  td  and  tdd  unknown;  f  -  all 
parameters  are  known. 
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parameters.  An  important  issue  of  practical  consequence,  which  remains  to  be  addressed  in 
future  work  is  the  effect  of  detector  response  time  on  parameter  estimation  of  echo  signals. 

The  double  point  target  provides  us  with  a  wide  parameter  space.  The  extension  of 
this  work  could  be  useful  in  introducing  the  natural  metric  in  parameter  space  for  the  proper 
evaluation  of  the  information  content  in  echo  signals.  The  bound  on  estimator  performance 
was  obtained  for  a  single  laser  pulse.  In  many  situations  of  practical  interest,  the  values  of 
Nsi,  NS2  and  ¥  will  fluctuate  from  pulse  to  pulse.  Evaluating  the  performance  of  a 
multipulse  system  is  certainly  of  great  interest.  A  few  minor  gaps  remain  to  be  filled  in  order 
to  make  this  analysis  complete.  This  includes  analysis  of  potential  benefits  of  pulse  wavefront 
design. 


One  of  the  results  obtained  was  that  for  large  relative  Doppler  shift  a  few  tens  of 
photons  are  required  to  measure  the  Doppler  shift  with  moderate  accuracy  (say  Avd  ~1).  To 
put  such  predictions  on  a  firm  footing,  one  has  to  construct  estimators  with  variances 
approaching  the  CRB.  It  is  easy  to  find  such  an  example  applicable  to  problems  in  which  fine 
parameter  estimation  is  of  interest,  e.g.,  alignment  and  calibration  when  the  uncertainty  of  the 
parameter  value  is  small.  The  form  of  locally  optimal  estimators  is  given  by: 
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This  estimator  is  for  small  displacement  of  A  from  A^.  It  is  a  biased  estimator,  however,  its 
bias  is  negligibly  small  for  small  displacements  from  Aq  and  its  variance  attains  the  CRB  at 

Ao- 


Another  important  class  of  estimators  useful  for  a  larger  class  of  applications  are  the 
maximum  likelihood  estimators  (MLE).  These  estimators  are  applicable  to  larger  parameter 
variations,  and  approach  the  CRB  for  large  signal  values.  The  methods  of  implementation 
for  the  MLE  estimators,  together  with  other  practical  considerations  and  types  of  estimators 
are  important  problems  which  should  be  addressed  in  future  work. 
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ABSTRACT 

A  novel  architecture  for  a  ladar  transceiver  is  proposed  to  achieve  high-resolution  imaging  and  precision  tracking  of 
complex  targets.  This  architecture  combines  a  conventional  imaging  receiver  (telescope)  with  a  noncoherent  sparse  array 
of  photon  buckets  into  one  system.  The  diffraction  limited  resolution  of  the  array  exceeds  that  of  the  telescope.  This  ladar 
architecture  is  designed  to  use  a  new  technique,  CRITIC,  which  is  based  on  combining  the  relevant  information  obtained 
simultaneously  with  the  telescope  and  the  array.  In  parucular.  the  CRITIC  technique  allows  one  to  overcome  several 
drawbacks  of  conventional  monostatic  ladars  and  opens  new  prospects  for  both  the  ground-  and  space-based  ladan.* 

1.  INTRODUCTION 

Fine  resoluuon  imaging  of  space  objects  and  precision  measurement  of  their  angular  position  is  of  interest  in 
numerous  suuauons.  Acove  laser  sensor  systems  are  considered  for  these  applications. 


Convenuonal  tracking  ladars  use  telescopes  (coherent  antennas)  to  image  a  laser  illuminated  target  onto  an  array  of 
photodetectors.  The  target  image  which  is  formed  using  a  single  laser  pulse  is  used  for  estimating  the  target  position. 
Typically,  the  image  resolution  and  angular  oosition  estimation  achievable  with  these  systems  are  on  the  order  of  the 
diffraction-limited  resolution  of  the  telescope1-2.  These  ladars  operate  at  the  short  wavelengths:  in  the  UV,  the  visible,  or 
SWIR,  in  order  to  take  advantage  of  the  diffraction-limited  performance. 


for 

for 


Simplicity  of  the  monopuise  ladars  is  their  main  advantage  for  space  applications:  direct  detection  method  is  used 
photosensing,  high  laser  spectral  purity  is  not  required1  -2,  and  the  staie-of-ihe-an  lasers  and  focal  planes  may  be  used 
illumination  and  detection. 


Conventional  monopulse  ladars  have,  however,  two  important  drawbacks:  (1)  They  are  not  robust  in  the  presence  of 
strong  specular  scattering  components  produced  by  the  target  glims  or  by  optical  countermeasures  such  as  comer  cubes^. 
(2)  The  telescope  images  may  be  highly  degraded  by  the  atmospheric  turbulence,  which  severely  limits  their  applicability 
as  ground  based  ladars  .  To  overcome  comer  cubes  as  the  countermeasures,  one  remedy  is  to  separate  the  transmitter  and 
the  receiver,  so  that  the  receiver  is  kept  away  from  highly  directional  comer  cube  reflections.  In  many  situations  such  a 
bistauc  configuration  is  r.ot  practical.  Also,  more  progress  is  needed  before  the  active  methods  for  controlling  the  figure 
of  large  telescopes,  and  the  adaptive  optics  designed  for  real-ume  compensation  of  atmospherically  induced  wavefront 
aberrations  will  meet  some  of  the  more  challenging  performance  requirements. 


Nonconvenuonal  active  imaging  sensor  systems  were  proposed  as  an  alternative  to  very  large,  expensive  imaging 
telescopes  using  adaptive  opucs.  We  refer  here  to  an  imaging  sensor  as  a  nonconvenuonal  when:  (1)  the  receiver  antenna 
is  noncoherent  and  consists  of  an  array  of  photon  buckets,  and  (2),  the  target  images  are  obtained  by  processing  the  target 
interferograms^. 
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The  reconstructed  images  approach  the  diffract  ion- limited  resolution  of  the  noncoherent  antenna  when  the  target  is 
within  an  isoplanauc  patch  of  the  atmospheric  turbulence.  Performance-degrading  laser  speckles,  are  the  well  known  side- 
effect  of  coherent  illumination.  They  can  be  reduced  by  ensemble  averaging  over  a  number  of  independent  realizations  of 
target  interftrograms  prior  to  the  reconstruction,  or  by  ensemble  averaging  of  images  reconstructed  from  the  individual 
interferograms. 

Nonconventional  imaging  systems  represent  an  attractive  alternative  to  conventional  imaging  systems  when  it  comes 
to  high-resolution  imaging.  Another  attribute  of  these  type  of  sensors  is  their  relative  robustness  to  the  presence  of  the 
strong  specular  components.  These  sensors  cannot  be  used  for  the  target  position  estimation,  however,  since  a 
noncoherent  array  of  photonsensors  cannot  measure  the  absolute  value  of  the  angle-of-amval. 

The  CRITIC  (Coherent  Radiauon  Imaging  and  Tracking  via  Nonuniform  Illumination)  technique,  which  we  propose 
in  this  paper,  makes  use  of  the  attributes  of  both  the  conventional  and  the  nonconventional  imaging  systems  and  is  free 
from  their  drawback*.  The  technique  is  introduced  in  the  next  section.  In  Section  3  we  provide  qualitative  examples  of 
the  CRITIC’S  performance.  The  last  section  is  a  summary. 

2.  THE  CRITIC 


In  this  section  we  describe  the  CRITIC'S  transceiver  and  its  operational  pnnciples. 
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Figure  1.  Schematic  of  the  CRITIC  Transceiver. 
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As  seen  in  Fig.  1.  the  system  receiver  consists  of  two  different  imaging  sensors:  one  is  a  conventional  imaging 
system,  in  which  telescope  images  are  sensed  at  the  focal  plane  with  an  array  of  photodetectors  and  another  is  a 
nonconventional  imaging  system  consisting  of  a  sparse  array  of  photon  buckets  for  sensing  the  target  interferogpuns  in  the 
aperture  plane.  The  diffraction-limited  resolution  of  the  telescope  is  typically  much  smaller  than  that  of  the  noncoherent 
antenna. 

A  single  transmitter  is  used  by  both  imaging  receivers.  The  transmitter  delivers  to  the  target  a  laser  pulse  with  a 
known  and  nonuniform  spatial  intensity  distribution  in  the  target  plane.  As  discussed  below,  the  known  beam 
nonuniformity  is  necessary  for  unambigous  reconstruction  of  target  images  from  the  interferograms  recorded  in  the 
aperture  plane.  The  transmitter  and  the  receiver  telescope  are  optically  aligned  with  high  precision.  This  alignment  is 
required  in  order  to  measure  the  beam  position  across  the  target  from  the  telescope  images. 

The  telescope  images  and  interferograms  are  collected  and  processed  in  parallel  for  each  illuminating  pulse. 
Whereas  the  telescope  images  are  produced  in  a  single  pulse,  the  reconstrucuon  of  images  from  the  interferograms  may 
require  more  than  one  beam  intensity  distnbuuon  across  the  target  and  therefore  more  than  one  pulse. 

The  mam  idea  behind  the  CRITIC  is:  (1)  use  the  telescope  images  for  measuring  the  angular  position  of  specular 
elements.  (2)  use  the  interferograms  to  produce  high  resolution  images  of  the  diffuse  elements  by  employing  a  technique 
recently  invented  by  us  called  self-reference  ambiguity-free  holography^,  (3)  derive  the  angular  positions  of  the  target 
diffuse  elements  relative  to  the  receiver  LOS  telescope  by  analyzing  the  images  obtained  with  both  receivers. 

2.1  TARGET  MODEL 


'Hie  complex  reflectance  of  the  target  can  be  modeled  as: 
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u> 


s  is  the  target  support,  a(r^)  and  A(r^)  arc  the  diffuse  and  specular  components,  respectively.  The  specular  component  is 
assumed  to  be  deterministic  and  is  due  to  the  target  glints  or  comer  cubes.  The  diffuse  component  is  due  to  an  optically 
rough  surface  and  is  modeled  as  a  complex  circular  Gaussian  process: 


<a(ro)>  •  0  (la) 

<»(r  )»(r  )>  *  0 
o  o 

<»(?  )«*(r  )>  ■  <la(r  )|2>  «<r  -  r  )  , 
o  o  o  o  o 

where  <  ■  >  stands  for  ensemble  averaging.  (In  Fig,  l  we  denote  the  diffuse  components  as  a)  and  aj.  and  the  specular 
components  as  A  j  and  A2  ) 

2.2  CONVENTIONAL  IMAGING  SYSTEM 

The  image  intensity  distribution  of  a  coherently  illuminated  target  in  the  focal  plane  rj  is  given  by: 


I(r  )  •  [dr  [dr  •  (r  )  a* (c  )  B(r  )B*(r  )  h(r  -  r  )  b* (r  -  r) , 

l  jOjoTolo  o  0  o  i  o  i 


(2) 


where  Bfr^  is  the  amplitude  of  the  incident  nonumiform  illumination  beam,  and  h(r£  -71)  is  the  coherent  point-spread 
function. 
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Consider,  for  example,  the  case  in  which  the  specular  components  are  seen  through  the  telescope  as  discrete  point-like 
sources,  and  are  several  orders  of  magnitude  brighter  than  the  diffuse  parts  of  the  target,  i.e.. 


•_<r  > 
T  o 


"  Z- 

n«l 


A  5  (r 
n  o 


(3) 


In  Refs.  1  and  2.  we  analyzed  methods  of  esumating  the  target  position  from  its  image  under  the  following  conditions:  the 
image  intensity  has  the  form  of  Eq.  2:  the  intensity  is  sensed  by  an  array  of  photodetectors  with  apertures  smaller  than  the 
Airy  disc  of  the  imaging  telescope:  and  the  accuracy  of  measurement  of  any  point  source  is  not  affected  by  the  presence  of 
others.  Under  these  condiuons,  we  have  shown  that  the  accuracy  of  measuring  the  angular  posiuon  Aj  of  a  point  source 
(see  Fig.  1)  is: 

A,  •  7<X/D,)MiT  ,  (4) 

*  2  g 

where  X  is  the  illumination  wavelength.  D2  is  the  diameter  of  the  telescope  aperture.  7  is  a  constant  which  depends  on  the 
point-spread  function,  and  Ng  is  the  expected  number  of  photocounts  from  a  glint  collected  by  all  photodetectors.  Eq.  4  is 
valid  in  the  limit  when  the  signal -generated  shot  noise  is  the  dominant  source  of  noise.  Typically,  practical  considerations 
such  as  quality  of  the  optics,  finite  extent  of  the  "point"  target,  and  nonuniformity  of  the  focal  plane  do  not  allow  one  to 
improve  the  angular  accuracy  bn  to  less  than  (XA^VOO  -  100)  merely  by  increasing  the  signal. 

In  the  CRITIC,  the  telescope  images  are  used  to  measure  the  angular  position  of  glints  rj,  and  their  strength  I  An  1^. 
The  angular  position  of  glints  is  measured  with  accuracy  better  than  the  diffracuon-limued  resolution  of  the  telescope 
(X/D2).  According  to  Eq.  4.  accuracy  100  umes  better  than  the  diffracuon-limued  resolution  can  be  achieved  when  a 
point-glint  produces  10.000  photocounts  (the  constant  7,  typically  does  not  exceed  unity)^.  According  to  Eq.  2, 
measurement  of  the  glint  strength  is  possible  when  the  illumination  beam  intensity  is  known  across  the  target.  In  the 
CRITIC,  the  knowledge  of  the  beam  intensity  distribution  is  achieved  through  optical  alignment  of  the  imaging  telescope 
with  the  transmitter  and  the  high  beam  quality.  The  accuracy  of  the  optical  alignment  should  be  at  least  as  good  as  the 
accuracy  of  the  angular  position  estimation  of  the  poim-glmts. 

2.3  NONCONVENTIONAL  IMAGING  SYSTEM 

The  CRITIC'S  nonconvenuonai  imaging  system,  shown  schematically  in  Fig.  1.  is  a  noncoherent  sparse  anray  of 
photon  buckets,  which  is  used  to  detect  interferograms  created  by  the  interference  of  all  target  components  under  the 
nonuniform  illumination. 


We  view  an  interferogram  of  a  complex  target,  as  a  hologram,  in  which  the  glints  serve  the  role  of  holographic 
references.  The  classical  off-axis  holographic  techniques  arc  not  useful,  however,  to  extract  the  images  from  such  an 
interferogram  since  the  reference  is  on  the  target  and  there  may  be  more  than  one  reference.  Under  such  condiuons  the 
conventional  techniques  will  produce  multiple  overlapping  images  with  unknown  angular  position  relative  to  the 
receiver^. 

The  self-reference  ambiguity-free  holography  which  wc  recently  proposed^  is  capable  of  extracting  the  images  of 
the  diffuse  parts  of  the  target  and  measure  their  angular  position  relative  to  the  nonuniform  illumination  beam. 

2.3.1  SELF-REFERENCE  AMBIGUITY-FREE  HOLOGRAPHY 
The  principle  of  the  CRITIC  holographic  technique  is  illustrated  in  Fig.  2  below: 
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Figure  2.  Schematic  of  Self-Reference  Holographic  Imaging. 

The  intensity  distribution  in  the  receiver  plane  m  a  distance  z  from  the  target- may  be  written 


Z  (r  )  a 


Ho  K 


,2»  - 

a  (r  )  *  (r  )  8(r  )  8  (e  )  axp  -ir-—  (*  -  *  )•  * 
o  T  o  T  o  o  o  [X*  oo 


where  B(r^)  is  the  amplitude  of  the  incident  nonuniform  beam,  and 

i_(r  )  a  a  (r  )  axp  [iii*  /  Xsl  .  (6) 

To  To  I  o 

Taking  the  Fourier  transform  of  the  interferogram  in  Ed-  (3)  gives 

r(r)  a  Ida  a?tr  )  a^Ie  -  r)  »(r)  B  (*  -  r) .  (7) 

In  the  above  equation  the  effect  of  finite  size  of  the  receiver  was  neglected.  The  effect  of  the  diffraction-limited 
resolution  on  the  reconstructed  images  was  investigated  qualitatively  and  is  illustrated  in  Section  3. 

Further  processing  of  information  contained  in  Eq.  (8)  depends  on  the  nature  of  target  reflectivity.  Consider  for 
simplicity  pointlike  glint  references.  In  this  case  the  target  reflectance  may  be  whiten: 

tt  _ 

•_(r>  •  S(r)  ♦  )>  I  »<?  -  5  >  ,  (•) 

T  *  -  ■  ■  n  q 

n«X 

and  the  Fourier  transform  of  the  interferogram  becomes 
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where 
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describes  interference  between  the  glints; 
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(11) 


is  the  interference  between  the  glints  and  diffuse  reflections;  and  lastly 


r 

oo 
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(12) 


is  the  autocorrelation  function  r.i  the  product  of  target  reflectance  and  if  e  beam  amplitude  distribution  in  the  target  plane. 

For  sufficiently  strong  glints  ( I  An  1  »  I  a(r)  I ).  the  term  FooVf)  may  be  neglected.  Since  glint  positions  and 
amplitudes  can  then  be  determined  with  high  precision  by  the  convent-onal  part  of  the  system,  the  term  Fgg  (f)  can  be 
removed  from  Eq.  (5)  by  an  appropriate  processing.  Therefore,  we  ,.vcd  to  consider  only  the  term  Fg0  (r). 

In  order  to  obtain  reconstructed  images  of  diffuse  targets  s  useful  to  consider  the  correlation  function; 

C  (r)  s  <f  (2)  ?*  U)>  .  (13) 

go  go 

Using  Eqs.  (la)  and  (11)  we  obtain 
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(14) 


This  is  an  equation  for  2N  unknowns  <  I  a(rj,  ♦?)  I  and  <  I  a(Tn  -7)  I  ^>.  The  unknowns  are  solved  by  generating  2N 
equations,  each  resulting  from  a  different  beam  position.^ 


Solution  of  He.  H  yields  2N  images  of  the  diffuse  target;  it  is  important  to  note  that  the  position  and  orientation  of 
each  image  is  known  w  >th  respect  to  positions  of  the  glints.  These  images  can  be  noncoherently  combined  to  enhance  the 
signal'io-noiJC  ratio. 
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13.2  DESIGN  CONSIDERATIONS 

Target  interferograms  are  spatially  band-limited  functions  with  the  highest  frequency  determined  by  the  target  size. 
The  highest  spatial  frequency  dictates  the  smallest  separation  among  the  photon  buckets  in  the  anay.  The  largest 
separation  defines  the  receiver  aperture  and  determines  the  diffraction-limited  resolution  of  the  array. 

In  the  CRITIC  architecture,  the  di  f  frre  tion  - 1  i  m  i  ted  resolution  of  the  array  is  one  or  two  orders-of-magnitude  larger 
than  that  of  the  imaging  telescope,  i.e.: 

(X/D3)  -  »<X/D2>  ,  (IS) 

where  D3  is  the  array  aperture  size,  and  m  is  always  larger  than  unity  and  may  be  as  targe  as  a  hundred.  The  resolution  of 
the  array.  A3,  depends  on  both  the  diffraction-limited  resolution  and  the  signol-to-noise  considerations. 

The  resolution  of  the  interferogram  depends  on  the  ability  to  measure  the  interference  fringes  between  the  diffuse 
and  specular  components.  The  signal  is  proporuonal  to  the  product  of  their  reflectivities  weighted  by  the  illumination 
intensity  across  them.  Assuming  orders-of-magnitude  larger  reflectivity  for  the  glint,  one  deduces  that  in  the  shot  noise 
limit,  the  glint  is  the  mam  source  of  the  noise,  and  the  signal  to  noise  rauo  is  proportional  to  the  square  root  of  the  number 
of  photons  from  the  diffuse  components  of  the  target. 

It  is  reasonable  to  expect  that  the  resolution  of  the  noncoherent  array  A3  can  not  exceed. 

A3  •  (X/D3)/(iT  ,  (H) 

where  Ns  is  the  total  expected  number  of  photons  from  the  diffuse  elements.  The  overall  performance  is  limited  by  the 
poorer  of  the  two  resolutions  t ^  and  A3. 

In  the  CRITIC  technique,  the  transmitter  is  used  to  code  the  target  reflectance  with  the  laser  illumination.  By 
illumtnauon  coding  we  understand  here  that  the  nonuniform  distribution  of  laser  intensity  is  known  across  the  target.  As 
we  stressed  in  the  section  2.3.1,  the  intensity  coding  is  the  key  for  unambtgous  reconstruction  of  the  target  images  from 
their  interferograms. 

The  longitudinal  coherence  of  the  laser  illuminator  must  be  at  least  twice  the  maximum  path  length  difference 
between  the  target  elements.  The  necessary  laser  power  can  be  deduced  from  the  system  resolution  requirement.  The 
sufficient  laser  power  is  larger  and  depends  on  specific  target  and  system  characteristics,  which  must  be  taken  into  account 
in  a  quantitative  analysis  of  the  sensor  system  performance. 

3.  CRITIC  PERFORMANCE:  AN  EXAMPLE 

In  this  section  we  present  results  of  computer  simulations  of  the  CRITIC  performance  in  the  strong  signal  limit, 
where  the  main  source  of  error  is  due  to  the  stochastic  nuiure  of  the  backscauered  field  arising  from  the  random,  target- 
induced  phase  modulation  which  results  in  the  well-known  phenomenon  of  the  laser  speckle. 

We  considered  a  model  of  a  complex  target  shown  in  Fig.  3A  illuminated  by  a  Gaussian  beam  shown  in  Fig.  3B. 
with  half-width  equal  to  the  target  length.  The  locations  of  "glints"  are  indicated  by  dots.  The  amplitude  of  the  'glint*  in 
the  center  is  10^  stronger  than  that  from  any  diffuse  element.  The  other  "glim"  is  10*  stronger.  The  weaker  'glint* 
models  a  strong  specular  component  from,  for  example,  a  satellite,  whereas  the  stronger  "glint"  models  a  comer  cube. 
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Figure  3.  Models  of  a  Target  and  Illuminating  Beam. 

A  CRITIC  transceiver  is  used  for  imaging  the  target  and  estimating  its  position.  Figure  4  shows  images  of  the  target 
obtained  with  the  conventional  telescope  receiver  considered  for  CRITIC  in  this  example.  Fig.  4 A  shows  the  target  image 
under  incoherent  illumination  and  in  the  absence  of  glints,  to  give  an  example  of  the  telescope  diffraction-limited 
resolution.  Fig.  4B  is  a  realization  of  the  telescope  image  when  illumination  is  coherent.  The  diffraction-limited 
resolution  is  represented  by  the  laser  speckle  size.  We  see  that  the  resolution  of  the  selected  telescope  is  too  poor  to 
recognize  the  target  from  a  single  realization  of  the  image,  due  to  the  noise-like  laser  speckle.  Figures  4C  and  4D  are  the 
telescope  images  of  the  target  in  the  presence  of  the  weaker  and  stronger  glints,  respectively.  As  expected,  the  telescope 
images  show  a  point-like  targets  whose  position  can  be  found  very  accurately,  but  not  a  trace  of  the  much  weaker  diffuse 
elements  of  the  target  can  be  found. 

However,  by  shifting  the  Gaussian  illumination  beam  to  a  different  position  and  using  the  two  different  resulting 
holograms,  we  obtain,  with  our  holographic  technique,  a  uirget  image  shown  in  Fig.  5C.  This  image  was  obtained  with  an 
uniformly  sampling  array  of  photon  buckets  with  an  aperture  four  times  larger  than  that  of  the  telescope.  Fig.  5D  shows 
the  reconstructed  target  image  after  averaging  over  100  realizations  of  the  laser  speckles  for  each  of  the  two  beam 
positions.  For  completeness.  Fig.  5B  illustrates  the  expected  failure  of  the  conventional  holographic  method  to  produce 
the  target  image.  The  importance  of  the  results  in  Fig.  3  is  that  they  show  target  images  which  are  produced 
simultaneously  by  the  two  sensors  employed  by  the  CRITIC  system.  Indeed,  one  can  find  the  glint  position  (Fig.  5A)  and 
have  high-resoluuon  images  (Figs.  SC  and  D). 

In  Fig.  6  we  show  similar  target  images  obtained  with  the  CRITIC  system  for  the  case  of  a  strong  glint  located  in  the 
middle  of  the  target.  The  noise  at  the  edges  of  Figs.  6C  and  D  is  an  artifact  due  to  the  tingle-precision  Fourier 
transformation  used  for  the  reconstructions,  while  the  dark  horizontal  lines  in  the  target  images  in  Figs.  6 C  and  D  are  the 
artifacts  of  the  specific  algorithm  used. 

In  Fig.  7A  we  show  a  single  realization  of  the  target  image  in  the  presence  of  two  glints:  the  weaker  at  the  top  of  the 
target  and  the  stronger  in  the  middle.  Images  in  Figs.  3A  and  6A  and  7A  were  simulated  for  the  same  illumination  and 
receiver  telescope.  Our  holographic  technique,  used  with  the  same  array  of  photon  buckets  as  before,  produces  images 
shown  in  Figs.  7B.  C  and  D  after  1. 10,  and  100  realizations  for  each  of  the  four  different  beam  positions  across  the  target, 
respectively. 

We  have  identified  several  promising  candidate  algorithms,  in  addition  to  those  used  to  produce  the  results  shown 
above.  We  continue  our  research  to  explore  the  CRITIC  technique  to  its  fullest  potential. 


# 


246  /  $Pf£  Vol.  1111  AcquiMion.  Tracking,  and  Pointing  HI  (1 969 ) 


A:  Non-coherent 
Illumination 


C:  Coherent 
illumination, 
Glint  s  102 
(1  realization) 


o  . 

4*. 


B:  Coherent 
Illumination, 
No  Glints 
(1  realization) 


0:  Coherent 
Illumination, 
Glint  s  IQ3 
(1  realization) 


Figure  i.  Telescope  Images. 
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B:  Conventional 
Holographic 
Reconstruction 
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Figure  5.  CRITIC  Images  for  Glint  Strength  »  10^. 
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Figure  6.  CRITIC  linages  lor  Gtini  Strength  *  10^. 
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Figure  7.  CRITIC  Images  for  Two  Glints  (10^  and  10^). 
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4  SUMMARY 


A  novel  ladar  system  is  proposed  to  achieve  high-rcsolution  imaging  and  precision  tracking  of  complex  targets  from 
a  ground-based  or  space-based  platforms.  The  system  is  based  on  the  CRITIC  technique  in  which  telescope  images  of 
glints  are  produced  simultaneously  wiih  holographic  images  of  the  diffuse  target  components  obtained  through  a  novel 
self-reference  holographic  imaging  technique  reported  by  us  recently.  In  the  technique,  the  glints  are  used  as  holographic 
references  for  the  much  weaker  diffuse  elements,  which  are  uniquely  coded  by  the  nonuniform  coherent  laser  illumination. 
By  integrating  into  one  system  the  imaging  telescope  of  modest  diffraction-limited  resolution  with  a  noncoherent  array  of 
photon  buckets  with  an  order-of-magnitude  larger  diffraction-limited  resolution,  the  CRITIC  overcomes  several 
drawbacks  of  the  conventional  monostatic  laser  radars  and  the  aperture  plane  interferometric  receivers,  while  preserving 
all  their  attributes. 
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